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TECHNICAL MEMORANDUM 13 83 

FINITE SPAN WINDS IN COMPRESSIBLE FLOW* 
By E. A. Krasilshchikova 


This work is devoted to the study of the perturbations of an airstream 
by the motion of a slender wing at supersonic speeds . 

A survey of the work related to the theory of the compressible flow 
around slender bodies was given in reference 14 by F. I. Frankl and 
E. A. Karpovich. 

The first works in this direction were those of L. Prandtl (ref. V) 
and J. Ackeret (ref. 23) in which the simple problem of the steady motion 
of an infinite span wing was studied. Borbely (ref. 25) considered the 
two-dimensional problem of the harmonically-oscillating nondeformable 
wing in supersonic flow by using integrals of special types for solutions. 

Schlichting (ref. 24) considered the particular problem of the flow 
over two-dimensional rectangular and trapezoidal wings. To solve this 
problem, he applied Prandtl' s method of the acceleration potential which 
he looked for in the form of a potential of a double layer. However, as 
shown later, Schlichting made an error and arrived at an incorrect result. 

In 194-3, Busemann (ref. 26) proposed the method of solving the prob- 
lem of the conical flow over a body by starting from the homogeneous 
solution of the wave equation. This method was modified by M. I. Gurevich 
who, in references 11 and 12, solved a series of problems for arrow-shaped 
and triangular wings when the flow, perturbed by the wing motion, is 
conical. The work of E. A. Karpovich and F. I. Frankl (ref. 13) was 
devoted entirely to the problem of the suction forces of arrow-shaped 
wings . 

In 1942, at a hydrodynamics seminar in Moscow University, Prof. L. I. 
Sedov proposed the problem of the supersonic flow over slender wings of 
finite span of arbitrary plan form. 

In response to this proposal of L. I. Sedov, there appeared in 
1946-47 a series of works by Soviet authors on the question of the super- 
sonic flow over wings of finite span. 

The first work in this direction was our candidate's dissertation 
(ref. 5 ) ) in which we found the effective solution for a limited class 


*Scientific Records of the Moscow State University, Vol. 154, 
Mechanics No. 4, 1951> PP* 1 81-239* 

The appendix represents a condensation made by the translator from a 
document "Modern Problems of Mechanics," Govt. Pub. House of Tech. Theor. 
Literature, (Moscow, Leningrad) 1952, pp. 94-112. 
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of harmonically-oscillating wings. In reference 6 we solved the problem 
for wing influences by "tip effect." Later works (refs. 15, 16, and 17 ) 
were devoted to the same problem. 

In reference 6 , using an idea of L. I. Sedov as a basis, we reduced 
the problem of the influence of the tip effect on harmonically-oscillating 
wings to an integral equation. 

The question of the flow over wings of finite span remained open for 
some time. 

At the start of 1947, there appeared works in which different methods 
were proposed for solving the tip effect problem which would be applicable 
to any particular wing plan forms. In reference 18 , M. D. Khaskind and 
S. V. Falkovich solved the problem, in the form of a series of special 
functions, for a harmonically oscillating triangular wing. Later, 

M. I. Gurevich generalized this method (ref. 19)- In reference 20, 

L. A. Galin reduced the problem of determining the velocity potential of 
an oscillating wing to the problem of finding the steady-motion velocity 
potential and gave a solution, in series, for the velocity potential of 
a rectangular, oscillating wing cambered in the direction of the oncoming 
stream. 

The methods, proposed by different authors, for solving the problem 
of the flow over wings of finite span do not permit the solution of the 
problem for any finite-span wing and may only be applied to a limited 
class of wings. 

Parallel developments in this direction were made by the foreign 
authors Puckett (ref. 2l) and Von Karman (ref. 22) who solved the problem 
of the steady flow over finite-span, symmetrical 1 wings at zero angle of 
attack. As is known, such wings produce no "tip effect" and the study 
of the perturbation of the airstream by their motion presents no mathe- 
matical difficulties. 

In references 6 , 7, and 8 we proposed a method of solving the finite- 
span wing problem by constructing and solving an integral equation which 
considered the wing plan form in both steady motion and oscillating 
harmonically. In reference 9 we generalized the problem to more general 
forms of unsteady wing motion by the method of retarded source potentials . 

Introducing characteristic coordinates we solved the integral equa- 
tion for wings of arbitrary plan form and represented the solution for 
steady wing -mot ion in quadratures and for the harmonically-oscillating 
wing in a power series of the parameter defining the oscillation frequency. 

The present work is a detailed explanation and further development 
of our papers (refs. 6 to 9) which were published in the Doklady, Akad. 
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Nauk, USSR. In this work we propose an effective method of solving 
aerodynamic problems of slender wings in supersonic flow. 

All the results and problems explained in this paper were reported 
by the author in 1947-48 to the USSR Mechanics Institute, V. A. Steklov 
Mathematics Institute, Moscow University, etc. 

In the first part of the work we find a class of solutions of the 
wave equation, starting from which we obtain the solution to the problem 
of determining the velocity potential of some wing plan form in unsteady 
deforming motion. The obtained solution contains the solution of the 
two-dimensional problem as a special case. In the same part of the work, 
we solve in quadratures the problem of steady supersonic flow over a 
wing of arbitrary surface and plan form. The effective solution for 
wings of small span is similarly given. We obtain formulas determining 
the pressure on the wing surface in the form of contour integrals and 
integrals over the wing surface. 

The author thanks L. I. Sedov for reading the manuscript. 


PART I 1 

1. SETTING UP THE PROBLEM 


1. Let us consider the motion of a thin slightly cambered wing at 
a small angle of attack. 

We will consider the basic motion of the wing to consist of an 
advancing, rectilinear motion at the constant supersonic speed u. Let 
be superposed on the basic motion, a small additional unsteady motion in 
which the wing surface may be deformed. 

Let us take the system of rectangular rectilinear coordinates Oxyz 
moving forward with the fundamental wing velocity u. The Ox-axis is 
directed opposite to the wing motion and we take the x,y-plane such that 
the z coordinates of points on the wing shall be small (figs. 1 and 2). 

We will consider the normal velocity component on both sides of the 
wing surface to be given by 

v n = A 0 + Ajjf jt + a] (l.l) 


^Results of Part I, sections 6 and 7 were found by the author in 
May, 19^7 at the Mathematics Institute, Akad. Nauk, USSR. 
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The first component defines the wing surface 

A 0 = -up 0 (1-2) 

where p 0 is the angle of attack of a wing element. The second compo- 
nent defines the additional unsteady motion of the wing. The functions 
Aq and Ap and a are considered given at each point of the wing 

surface . 


We will assume that the fluid motion is irrotational and that there 
are no external forces . 

The velocity potential of the perturbed stream cp(x,y,z,t) is 
represented in the form 

Cp(x,y,z,t) = cp 0 (x,y,z) + q>l(x,y,z,t) (l-3) 

where the potential cpQ corresponds to the basic steady motion of the 
wing and the potential corresponds to the additional unsteady motion. 

Thus the projections of the velocity v of the fluid particles on 
the moving Oxyz coordinates are determined by 


dcpQ dq>l 

dx dx 


d^o 

dy dy 



dcpQ dcpj 

dz dz 


The functions cp Q and cp^ and their derivatives will be considered 
first-order quantities and second-order quantities will be neglected. 

With these assumptions it is known that the potential cp^ satisfies the 
wave equation which in the moving axes is 


(a 2 - u 2 ) 


d 2 cpi 2 8 2 cp- 


+ a' 


1 + a 2 


dx 2 dy 2 

and the potential cpQ satisfies 

/ P P\ a2( P 0 2 

(a 2 - u 2 ) + a 2 

dx 2 


d 2( Pl _ ^1 _ 2u d2 

dz 2 dt 2 dtdx 


q>L _ 


= 0 


(1.4) 


d 2 (p 0 [ ^2 8 2 Tq _ ? 

dy 2 dz 2 


(1.5) 


where a is the speed of sound in the undisturbed stream. 

A vortex surface, called the vortex sheet, trails from the side of 
the wing surface opposite to its motion. Just as on the wing surface 
the velocity potential undergoes a jump discontinuity on this sheet. 
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We represent the projection of the vortex sheet on the x,y-plane as the 
semi-infinite strip £j_ (fig. l) extending along the x-axis to infinity 
from the trailing edge of the wing. 

Let us establish the boundary conditions which the functions cp^ 
and cpj_ satisfy. 

Let us transfer the boundary conditions on the wing surface parallel 
to the z-axis onto the projection £ of the wing on the x,y -plane, 
which is equivalent to neglecting second-order quantities in comparison 
with first-order ones. Therefore on the basis of equation (l.l) we obtain 
the streamline condition 


^ = A 0 (x,y), — = A x (x,y)f [t + <x(x,y)l 

oz oz 


( 1 . 6 ) 


which must be fulfilled on both the upper and lower sides of Z. 

The kinematic condition, which expresses the continuity of the normal 
velocity components of the fluid particles, must be fulfilled on the dis- 
continuous surface of the velocity potential and on the vortex sheet. 

We transfer the condition on the vortex sheet parallel to the z-axis 
onto its projection Zj_ on the x,y-plane which is again neglecting second- 

order quantities. Therefore we have the conditions 


8cp 0 


dcpo 


691 


dcpi 

dz 

z=+0 

dz 

N 

II 

1 

O 

V. 

dz 

z=+0 

dz 


( 1 . 7 ) 


to be fulfilled on Z^_. 

Furthermore, the dynamic condition which the potentials 9 q and 
cp-^ satisfy must be fulfilled on the vortex sheet. 

Since the pressure remains continuous on crossing from one side of 
the vortex sheet to the other, then from the Lagrange integral 


P = _ le _ u £2 - i itesef + + (*£\ 


dt 


dx 


dx 


\dy) 


K dzJ 


V 

> + r(t). 

hd ' 
11 

j 
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Keeping equation (1.3) in mind and neglecting second-order quantities, 
we obtain 


dcp 0 




dcpi + 

dcpi 

u 


a a + 

dcpi 

u 

dx 

z=+0 

dx 

z=-0 

dt 

dx 

z=fO 

dt 

dx 


which must also be fulfilled on Z]_. 

After boundary conditions ( 1 . 6) and (1.7) are established, we 
correctly consider that, to the same degree of approximation, the surface 
of discontinuity of the velocity potential - the vortex surface - lies 
entirely within the x,y-plane. Therefore, the functions cpQ and cp^ 

are odd functions in z 


= -<P 0 (x,y,z), qp x (x,y,-z,t) = -cp 1 (x,y,z,t) (1.9) 


Combining equations (l.8) and (l-9) we conclude that the functions 
9 ^ and cp^ satisfy the respective conditions 


dq^ 

8x 


= 0 , 



dcp-i 

u = 0 on Z]_ 

8x 


( 1 . 10 ) 


Since the motion of the wing is supersonic, the medium is disturbed 
only in the region bounded by the respective disturbance waves represent- 
able by a surface enveloping the characteristic cones with vertices at 
points of the wing contour. Ahead of this surface - in front of the wing - 
the medium is at rest, therefore, the velocity potential is a constant 
which we assume to be zero. Hence we have the condition on the disturb- 
ance wave 


cp 0 (x,y,z) = 0 , <Pi(x,y,z,t) = 0 (l.ll) 

The potentials cp Q and cp]_ are continuous functions everywhere 
outside the two dimensional region Z + Z;j_ and, as was established, are 
odd in z, therefore, in the whole x,y-plane outside of the region £ + Zq 
where the medium is perturbed, the following conditions are satisfied: 


cp Q (x,y,0) = 0, 


9q( x ,y>0,t) = 0 


( 1 . 12 ) 
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The region where equation (1.12) is satisfied is denoted in figure 1 
by s 2 and 2 2*‘ 

Thus the considered hydrodynamic problem is reduced to the following 
two boundary problems: 

I. To find the function cp-^x^^t) which satisfies equation (l.4) 
and boundary conditions (l.6), (l.lO), (l.ll), and (1.12). 

II. To find the function cpQ(x,y,z) which satisfies equation (1.5) 
and boundary conditions (l.6), (l.lO), (l.ll), and (1.12). 

Since the functions cpQ and cp^ are antisymmetric functions rela- 
tive to the z = 0 plane, it is sufficient to solve the problem for the 
upper half plane. From the solution of boundary problem I it is possible 
to obtain the solution of II if the function f in the first be considered 
a constant equal to unity, and Aq replaces Ap. 

2. VELOCITY POTENTIAL OF A MOVING SOURCE WITH VARIABLE INTENSITY 


1. Let us construct a solution of equation (1.4) as the retarded 
potential of a source moving in a straight line with the constant velocity 
u and having an intensity which varies with time according to fp(t) . 

Let us consider the infinite line along which, at each point from left 
to right, sources with velocity u start to function one after the other 
with the variable intensity q = f()(t - tp)fp(t). The law of variation 

of the function fQ is the same for all the sources if the initial 
moment of each source is considered to be the moment when it came into 

p 

being. 

The function fp has the same value for all the sources at each 
instant. Let a source at an arbitrary point of the O'x'-axis be acting 
at time tp (fig. 3) • The retarded potential of the velocity at the 

point M as a result of such a system of sources is represented in the 
fixed coordinates by 



2 Prandtl (ref. 3) considered an analogous problem with q = fo(t - tp) . 
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where A is a constant with the dimensions of a velocity. The limits 
of integration t-,' and t]/’ take into account those sources which 
affect M at time t. The origin of the fixed coordinates O' is 

placed at the point at which the source started at t = 0. 

Introducing the new variable of integration t = a(t - tp) - r and 

transforming to the coordinate system x = x' + ut, y=y', z = z ' 

which is moving forward in a straight line with the velocity u, we 

transform equation (2.1) into 



jf it is assumed that u > a then the velocity potential at M(x,y,z) 
is the sum of the expressions (2.2), with the minus sign in front of the 
radical taking into account the effect of the sources in the strip AC 
on M and with the plus sign taking into account the sources on CB. The 
smaller root of the radicand is taken as the upper limit of integration 
x it is easy to see that in this case both roots are real, positive 

quantities (fig. 3) • 


On the basis of expression (2.2) we now construct a velocity potential 
at M from the sources moving with speed u > a which have an intensity 
which varies with time as f x (t) . The derivation remains valid if the 
additive constant a x is added to the argument t of the function f x . 

Putting the sources at the origin, we find the velocity potential from ^ 
equation (2.2) by considering the interval of integration from 0 to t x 

to be vanishingly small. Then, neglecting the term and - putting 


A 

a 


JO 


f q dx = C where C is a constant, we obtain the desired solu- 

1 


tion for equation (1.4-) in the general form 
f^\t + - 


ux 


<P*(x,y,z,t) = C 


u2 - a 2 u2 


1x2 - 


ui 

Va2 


(y2 + z2) 




iVy 2 + z 2 



(2.3) 


MCA TM 1385 


9 


Let us note that each component of the arbitrary function fq as 
well as the constant C and aq in equation (2.3) is separately also 
a solution of equation (l.A) . 

In equation (2.3) putting ay = 0 and the velocity of motion of the 
source u = 0, we arrive at the well-known solution for a spherical wave. 

If the velocity of motion of the source is u < a then to obtain 
the retarded potential of a moving source the right side of equation (2.3) 
must be limited to the first component. 

Considering the function fq in equation (2.3) to be constant, we 
arrive at the Prandtl (ref. 3) solution for the retarded potential of a 
moving source of constant intensity 



2. It is possible to obtain, by the same method, the velocity 
potential of a source with the variable intensity fq(t) moving 
arbitrarily. 


For example, in the case of rectilinear motion of the source when 

aFqtt) 


the motion is given by X = Fq(t), Y = 0, Z = 0 and when 


dt 


> a. 


that is, the motion of the source is supersonic, the velocity potential 
of the source at the origin of a coordinate system moving with the source 
is 




Cfj(ti) 


y jx + Fi(t) - ^ + y2 + z2 _ jjc + F^(t) - F^(t)J 


dFxCtg) 

dtx 


Cfx(ti*) 


£ + F x (t) - F-^tx*)] 2 + y2 + z 2 _ j^x + Fj_(t) - F 1 (t 1 *)] 


dFx(ti*) 

dt x * 


( 2 . 1 ) 


where the parameters tq = tq(x,y,z,t) and tq* = tq*(x,y,z,t) are 
real roots of 


l(t - tq) - ^[x + Fq( t) - Fq(tq)J 2 + y 2 + Z 2 = 


0 


(2.5) 
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If 


[dFp(t)/dt] 


< a , i.e. 


the source velocity is subsonic, then to obtain 


the velocity potential one must be limited to the one component in equa- 
tion (2.1+) which corresponds to the smaller of the values of the parameters 
t-L and tp*. 


The function expressed by equation (2.4) 
tion with variable coefficients 


satisfies the linear equa- 


dF x (t) 

2 

8 2 a 2 a 2 

l + a 2 LfP + a 2 -1 - 

a 2 cp 

dt _ 

j 

1 8x 2 8y 2 8z 

at 2 


; "1^ a 9 

dt chcdt 


a' 2 Fi(t) \ 

dt 2 


(2.6) 


If the source moves with constant acceleration as F-i(t) = -ut - 

x 2 

(where b is a constant) then equation (2.5) is an algebraic equation 
of the fourth degree in t]_ with two real roots. 

Formula (2.4) contains the Lienard-Weigert (ref. 27) formula as a 
special case when the source intensity is constant. 


3. DERIVATION OF THE BASIC VELOCITY POTENTIAL FORMULA 

1. We apply a solution of the form (2.3) of the wave equation (l.4) 
to the above-mentioned boundary problem I . 

At each point of the x,y-plane let us place sources with the poten- 
tial cp*. Hence, we will consider C and in equation (2.3) functions 

of points of the x,y-plane and we will replace ap by a and fp by f. 

As a consequence of the linearity of equation (l.4), its solution 
is a function cpp expressed by 


<p 1 (x>y> z > t ) = J J 

s (x.y, z ) 


// «!,.) 
S(x,y,z) 


t + a(£,Ji) - U g X * ■ - 2 a 2 J(* - s ) 2 - * 2 (y - n ) 2 - k2 * 2 \ 

: U 2 - q2 u 2 - a2 V I dl)d | + 

^(x - s ) 2 - k 2 (y - I )) 2 - k 2 z 2 




- eP- - b?- 




f (x - i ) 2 - k 2 (y - t )) 2 - k 2 z 2 


di]d£ (5-1) 


\J(x - |) 2 - k 2 (y - t ]) 2 - k 2 z 2 



where k = 


1. 
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The region of integration S(x,y,z) is that part of the x,y-plane 
which lies within the characteristic fore-cone of equation (l.U) from 
the point with coordinates x,y,z (fig. k) . 

The solution of equation (3.1) will give the velocity potential 
arising from the additional motion of the wing if C(x,y) is determined 
from the boundary conditions of the problem on the x,y-plane. 

Let us introduce the new variable of integration 9 into equa- 
tion (3-l) in place of t) 


t) = y - ~ \J(x - I) 2 - k 2 z 2 cos 9 
Then equation (3*l) becomes 


(3-3) 


cp]_(x,y,z,t) = J J CU,y - i )j(x - i) 2 - lz‘ 

S(x,y,z) l 


2 z 2 cos 0 V x 


f -S t + a 


l,y - J ^|(x - |) 2 - k 2 z 2 cos 0 


u ( x ~ l) S 

2 2 2 2 
vr - a u - a 


\j(x - |) 2 - k 2 z 2 sin 9 >d0d| + 



z) 


C< 


t>y 


i \J(x - |) 2 - k 2 z 2 


cos 9 / x 


f J t + a 


t>y 


-i'F* - 5)1 


A 2 


COS 0 


u(x - 5) + a 

2 2 2 2 
vr - &r - or 



, 2 2 
k z 


sin 0 


i 


d9d| 


(3-^) 
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Let us note that for any point M(x,y,z) of space it is possible 
to isolate from the region S(x,y,z) a region S' in which the vari- 
able of integration has the limits 

x-kz=|=C', 0 ^ 8^jt 


or 

n 2 = y - £ |/u - 1) 2 - k 2 z 2 < ri < y + i y^(x - 1) 2 - k 2 z 2 = 

where C' is a constant satisfying the inequality C' < x - kz. In the 
remaining region S - S' the limits of integration either do not depend 
on z or depend on z only in the combination kz^. 

Differentiating equation (3.I+) with respect to z we find the rela- 
tion between c(x,y) and a(x,y) and the normal derivative of the 
velocity potential dcp^/dz at any point of the x,y-plane 


C(x,y) 


i <jf|t + a(x,y)| 


_dz 


z=0 


(3.5) 


Comparing equation ( 3 * 5 ) with equation (l.6) we conclude that on 
the wing 


C(x,y) = - Ap(x,y) ( 3 - 6 ) 

i.e., the function C(x,y) is given. 

Therefore, the velocity potential may be computed from equa- 

tion ( 3 . 1 ) by taking equation (3.6) into account for those points M(x,y,z) 
of space for which the region of integration S(x,y,z) does not extend 
beyond the limits of the wing. 

If the leading and trailing edges of the wing are given by x = \|r(y) 
and x = XpCy)/ respectively, and if, therefore, i|r and X 1 satisfy 
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g- i( yl 

ay 


< cot a* 


(3-7) 


d*i(y) 

ay 


< cot a* 


(3-8) 


(where a* is the semi-vertex angle of the characteristic cone) on the 
leading and trailing edges of the wing, respectively, then in particular, 
equation (3.1) yields the effective solution of the problem of finding 
the velocity potential cp-^ everywhere on the wing surface because in 

this case the region of integration S does not extend beyond the wing 
for any point M(x,y,0) on it (fig. 5) * 


Also, in particular, equation (3*l) gives a solution of the plane 
problem if C and a axe considered as functions of one variable - 
C = C(x) and a = a(x) - and the variables of integration in the region 
S are considered to vary between 


0 < | < x - kz 

tj 2 = y - i \J(x - |) 2 - k 2 z 2 < T ) < y + 1 \J(x - |) 2 - k 2 z 2 = 

(3.9) 

where and are as defined previously. 

Considering f in equation (3«l) a constant and taking into account 
equation (3.5)* we obtain the fundamental formula for the velocity poten- 
tial cpQ specified by the basic steady motion of the wing 


<P 0 (x f y,z) 


i r r 


it ,J J . 

S(x,y,z) 

dz 



d-Tjdg 

\/(x - 5) 2 - k 2 (y - T )) 2 - k 2 z 2 


( 3 . 10 ) 

Formula (3.IO) contains, as special cases, the results of Prandtl 
(ref. 3) y Ackeret (ref. 23 ), Schlichting (ref. A) when the wing surface 
is a plane and when the leading edge is a straight line perpendicular to 
the free stream. 
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4. HARMONIC OSCILLATIONS OF A WING 


1. Let us turn to the case when the additional motions of 
are harmonic oscillations, i.e., on the wing equation (l.6) is 


the wing 
given as 


d<J>L 

dz 


R.P. A;L(x,y) 



+ a 


(x,yj| 


= R.P. A P (x,y)e 


in>t 


(4.1) 


where A 2 (x,y) defines the amplitude and initial phase of the oscillations. 

Using the obvious relation e i0 + e“ i0 = 2 cos 0 and equation ( 3 . 5 ) > the 
basic formula for the velocity potential ( 3 . l) is represented as 


<P 1 (x,y,z,t) = - i e px 

x 31 


dz 


e“^cos 


A \j(x - £) 2 - k 2 (y - T ]) 2 - k 2 z 2 


S(x,y,z) z -0 - l ) 2 - k 2 (y - T ]) 2 - k 2 z 2 


dT)d| 


( 4 . 2 ) 


where 


and 



u 2 - a 2 

Keeping the second inequality of equation ( 3 . 9 ) in. mind, let us 
compute the inner integral after which we obtain a solution of the prob- 
lem for a wing of infinite span 


cp x (x,z,t) = - £ 




-oc-kz 


dz 


,-PI- 


- 1) 2 - 


v 2 2 

k z 


46 


z=0 


(4.3) 

where Iq is the Bessel function of zero order. 

By means of equation (^. 3 ) the velocity potential may be computed at 
those points of the x,z-plane for which the interval of integration on 
the Ox-axis does not extend beyond the wing, i.e., at those points of the 
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x,z -plane not affected by the vortices trailing from the wing because 


the function 


dq>i 

dz 


is given only on the wing. In order to compute the 


velocity potential at any point of the x,z-plane by equation (I4-.3) it is 

dcp-, 

necessary to determine — =, using equation (1.8), everywhere on the 

dz 

Ox-axis outside the wing. 


Let us express, by equation (^.3), the velocity potential cp^ for 

any points lying on the Ox-axis outside the wing, which, according to 
equation (l.8), equals on the Ox-axis everywhere outside the wing 


where 


^(xjt) = R.p. 

v--*E 

u 


( 4 . 4 ) 


and l is the abscissa of the trailing edge. Then we obtain the integral 
equation 


'-^pi e" p ^I 0 |A(x - {)|d| = -kcf^e 


« ~- 3 x 


0 


e" PI Io 


A(x - |)[ d| 


! z=0 


(4.5) 

dcp 

which - — satisfies on the Ox-axis outside the wing. In reference 
oz 

we solved such an integral equation. The inversion of equation (I1-.5) is 



dF*(x) 

dx 



F^( 5 )It |A(x - O 


x - i 


( 4 . 6 ) 


where F* denotes the right side of equation ( 4 . 5 ), the known function, 
and where 1 ^ is the Bessel function of first order. 

Therefore, keeping equation ( 4 . 6 ) in mind, we can calculate the 
velocity potential at any point of the x,z-plane by equation (4-3). 
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The problem considered in this section was solved and explained 
in reference 5 from another point of view. 


5 . INFLUENCE OF THE TIP EFFECT 


1. To calculate the velocity potential according to equation (3-l) 
and also through equation (3*10) or (L.2) for those points M(x,y,z) of 
space for which the region of integration S extends outside the limits 
of the wing surface, it is necessary to determine the normal velocity 


component 


8q>L 

8z 


everywhere in the region of integration S 


from the 


boundary conditions of the problem on the z = 0 plane. 


Let us consider the case when the region of integration S inter- 
sects the wing surface and the region E^ lying outside the wing and 
outside the region of the vortex system from the wing. Region E^ 

(fig. 6 ) is part of the region E 2 defined above. That is, let us con- 
sider the case when the wing tips - the arcs ED and E'D' of the wing 
contour - act on the point M(x,y,z) or so to speak, the influence of 
the "tip effect" and not the influence of the vortex sheet trailing from 
the wing surface . 

The point E on the leading edge is defined so that condition (3-7) 
is fulfilled to its left and violated to its right. The point E' is 
similarly defined. The points D and D' are, respectively, the right- 
most and leftmost points on the wing contour as shown in figure 6 . 

/ ^ dcp l 

Let us construct the integral equation for C(x,y), connected to - — 

oz 

by relation ( 3 - 5 ) > in ^ 5 - 

Let us select the velocity potential at any point N(x,y,0) 

lying in E^ by means of equation ( 3 . 1 ), equal to zero everywhere in E^ 

according to equation (l.l2). The region of integration S(x,y,0) is 
divided into two parts, as shown in figure 7; the region s(x,y) is 
that part of the wing falling in the Mach fore-cone from N(x,y,0), and 
the region o(x,y) is that part of E^ lying in the same fore-cone. 

According to equation (3*6) C(x,y) is given in s. In o, c(x,y) is 
unknown. We therefore arrive at the integral equation which C(x,y) 
satisfies in E^. 
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IT 


cr(x,y) 


C(|,T])K(|,T];x,yjt)dT|d| = F(x,y,t) 


(5-1) 


where the kernel is 


K = U,n;x,y,t) = 


fit + a(|,ri) - ^ 4 - 2 a g l/(x - |) 2 - k 2 (y - r,) 2 j 

L - &r ur - 1 + 

\/(x - i ) 2 - k 2 (y - n ) 2 


•{» * - ife# * nrhg ^ - 5)2 ~ t2(y ~ " )2 } 




- t ) 2 - * 2 (y - 'll 2 


(5-2) 


and the known function 


F(x,y;t) = — .JJ ^ A 1 (t,n)Kit,wx.,y}t)dvte (5-3) 


Kx,y) 


If the characteristic coordinates are introduced 


= x 


Xq - k(y - yo) , y L = x - x Q + k(y - y 0 ) , z ] _ = kz 


(5-^) 

(where xq anil yQ may he any numbers) then integral equation (5*l) is 
simplified and in some cases this integral equation is easily inverted 
as will be shown below. 


6. SOLUTION OF THE INTEGRAL EQUATION FOR A HARMDNICALLY OSCILLATING WING 

1. If the additional motions of the wing are harmonic oscillations, 
i.e., the condition on the wing is given in the form of (U.l), then 
equation (5-l) becomes 





COS 


[a \/(x - |) 2 - k 2 (y - T])‘ 
\/(x - g ) 2 - k 2 (y - t )) 2 


— dT]d5 = F(x,y) 


( 6 . 1 ) 
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where the function 9(x,y) = 
function is 


dcp 


8z 


z=0 


S "P X in a and where the known 


K*,y> . - rr A(i,n) cos[a ft* - if - **(y - 1)11 a , as ( 6 . 2 ) 

, JJ , /(* - t) - k 2 (y - -l ) 2 

s(x,y) 


where A(x,y) = • 


Scp i 

8z 


-px 


z=0 


in s. In order to solve this integral 


equation we introduce the characteristic coordinates x^, y^, z^ with 
origin at "0" hy means of the formula 


X]_ = x - ky, yq = x + ky, Zq = kz 


(6.3) 


In the new coordinates the variables of integration in a will vary 
between the limits 


*E “ S 1 - *1' *(?l) ” \ _ y l 


(6.4) 


where yq = ty(xq) is the equation of the wing tip - the arc ED of 

the wing contour - in the transformed coordinates, and Xg is the 

abscissa of E defined in section 5 i n these same coordinates (fig. 8). 
Equation (6.1) is transformed to 

- 4 \«i - y i) 

(6.5) 


r 

X E 


y l 

0 -» (I 
t(il) lV 


1 ' 


\) 


3 [*lfo. - e i)( y i - ’’i 

/( x i - h)( y i - M 
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where the function 


e i( x i>yi) 


8cp x 



e 


P ( x i+yi) 

2 


and where the known function is 


F^yj.) - - ff Ai(£j..’» 3 .) C 0 S [V( x i - s jjhi - id] at)i iSi 

s(x a ,yi) it* 1 ' s h( y i - ’’i) 

w 


A l» 


8zi 


P ( x i+yi) 
2 


z i=° 


(6.6) 


8cp x 

Let us note that the normal velocity of the perturbed flow — 

O Zj- 

is related to 8cp^/8z]_ by 


52 = 

8z 8z]_ 

For brevity, the index "l" will be left off the independent variable 
everywhere from now on. 

2. Let us look for a solution of equation ( 6 . 5 ) in the form of the 
power series 

00 O 

0 (x,y;A) = ) G^Cx^) A n ( 6 . 7 ) 

n=0 


Into both sides of equation ( 6 . 5 ) let us introduce 


cos 


X\/(x - tXy - 1)1 = in (x - i) n (y - T|) n (6-8) 
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Keeping the absolute convergence of equations (6.7) and (6.8) in 
mind, we multiply them term by term with the result 

e(l,n)cos[A /(x - 4 ) ( y - n) ] 


= 11 A 211 !! ( '^- k [(X - 4 ) (y - n)] n - k 0 (4,11) (6.9) 

n=0 k=0 L 2 ( n “ k )J- L J ^ 

Substituting equations (6.7), (6.8) and (6.9) into equation (6.5) 
the latter becomes 


x y 


Xg t(l) 


^ T 2 ' j ' ~ 1 ' > k d, 0 2k (M)[(* -5) (y - l)] 

n=0 k=0 [2(n - k)Jl L J 


n-k-i 


dii d 4 


>(x,y) 


A(l,n) f: ■ 0 1 ) n+1 - s> (y - n)] 2 an 

SO (2n)l 


n- 


d| 


( 6 . 10 ) 


Taking into account the uniform convergence of the series in both 
sides of equation (6,10) with respect to the variables g and T] we 
integrate term by term 


oo k=n n — k x x 

r- , 2 a y- (-D n r r 

fco k=0 f2(n - k)1 : J J 

^ +(i) 


n-k-i 

e 2 k [( x ” (y - 2 dT1 d ^ 


n=0 


2n (-1) 


n + 1 


4 


(2a)! 


JJ A(4,ti)[(x - 4) (y - ri)] drj d4 (6.11) 


s(x,y) 
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In equation (6. 11) equating coefficients in identical powers of A 
we obtain the integral equation which the functions ®2n^ x -’ y ) satisfy 


p Py 

Ml) 


9 - (M) CTI=T - F n<^> 


where 


n-1 k 

F„(x,y) = f n (x,y) + ^ f„(x,y) 




n' 


(6.12) 


(6.13) 


where, in its turn, 

fn(x ’ y) =i "Sr~ ^ A(M)[(x ' 5)<y " ^ ^ 45 (6 - l4 > 

s(x,y) 

and 

V f 1 \ Pv Pv r- o’ 

fn<X,y) = [2(n - k)] ! L i(S) %k<M) ^ " l)(y " n)] 4,1 41 

E (6.15) 

Ir- 

from which the functions f n are defined for k 5 0 and n > 0. Let 
us note that the right side F n (x,y) of equation (6.12) depends, for 
0 2n , on the coefficients 0 2k hut only for k = 0,1,2, .. ,,n-l. There- 
fore, if we find 0 O , 0 2 , 0^, ..., e 2 ( n _i), then F n (x,y) is a known 

function in the equation which the coefficient 0 2n in the general term 

of series equation (6.7) satisfies. For n=0 the right side in equa- 
tion (6.12) 


Fq (x,y) = %(x,y) 



s(x,y) 


is a known function of x and y. 

Let us solve equation (6.12) for 


A(i,n) d,) dg - 

|/(x - S)(y - T) 

e 2n^ y h 


(6.16) 
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The two dimensional integral equation (6.12) is equivalent to the 
two homogeneous integral equations 

P X 0* (f y) 

j- dl = F n (x,y) ( 6 . 17 ) 

JXj, * x - 5 

and 

n y 

4 ( 1 ) e 2n(S,,l ) an = e^d.y) ( 6 . 18 ) 

vy - n 

each of which reduces to an Abel equation. 

Using the inversion formula of the Abel integral equation and 
observing that for any n functions F n (xg,y) = 0 hence the solution of 

equation ( 6 . 17 ) for the function 02 n (x,y) is 



r x F ne (!,y) 

^ / x - 6 


d£ 


(6.19) 


Let us turn to equation (6.l8). We denote the parameter £ by x, 
and again using the inversion formula for the Abel equation and keeping 
in mind that according to equation ( 6 . 19 ) the right side X, i(x)[ of 

equation ( 6 . 18 ) for y = i|r(x) is different from zero, the solution of 
equation ( 6 . 18 ) for 8 p n is 


e 2n (x ,y) 


1 ® 2 n [ x ?^( x )1 + 2 . py e * 2 nr ]( x > T )) 
31 fy ~ *U) * J *(x) j/y - n 


(6.20) 


Substituting in equation (6.20) in place of 0* n (x ;1 y) its value from 
equation ( 6 . 19 ) we obtain the solution of equation (6.12) in the following 


form: 


02n( x ^y) = 


it 2 i/y - y(x) Jx e /x - i 


d£ + 


Jx E J 1 | f ( x ) /(x - £) (y - rj) 


dq d£ 


( 6 . 21 ) 
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Thus, according to equation (6.21), we can evaluate successively, 
the coefficients 0 q, ®2k> etc * 

Formula (6.21) shows that all the coefficients (n=0,l,2, . . . ) for 
y = \|r(x) , i.e., on the wing tip ED, become infinite as r-1/2 

where R 

is the distance of the point (x,y) from ED. Therefore, the velocity 
of the perturbed stream becomes infinite as the specified order on the 
wing tips, approaching from outside the wing. 

It is possible to represent the inversion (6.21) of (6.12) as 


3 2n 


(x,y) = -i- 


8 2 


z 


F n (M) 


*2 dxdy J XE J t ( x ) - 0(y - n) 


dq dg (6.22) 


which can be confirmed without difficulty by direct differentiation with 
respect to the parameter. 

Therefore, the solutions of integral equation ( 6 . 5 ) are constructed 
in the form of the absolutely convergent series ( 6 . 7 ) for any value of 
the parameter "X . 

The coefficients 6pn(x,y) are expanded in the series 


0, (x,y;X) 0 2n^ x,y ^ ^ 

n=0 


(6.23) 


We find the function 0*(x,y) 



z=0 


e 


-0 


(x+y) 


2 


in. 


(fig. 6) lying off the wing to the left, from equations (6.21) or (6.22) 
by replacing in the latter the function y(x) by ^(x) (where 

y = ^(x) is the equation of the arc E'D 1 of the wing contour - 

the left wing tip) and interchange the role of the coordinates. 
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3- Let us consider a wing of small span. Let the characteristic 
cones from Ep and Ep‘ intersect the wing as shown in figure 9* The 

points E 1 and E]_' are defined just as are E and Ep in section 5 . 

Let us divide the x,y-plane where the medium is perturbed into the 
regions So, Sp, S 2 , • • • y S n , .... 

The region S n is the M-shaped region lying within the character- 
istic aft-cones from E n and E ^ 1 (or within one of them) and outside 
the characteristic aft-cones from E n ^.]_ and E n+ p' . In its turn, we 

divide the part of the x,y-plane lying to the right and left of the wing 
into the strips bp, (T 2 , • • *, G n , . . . and dp', G<^ , ' ‘ ' ’ 

a n ', . . ., respectively. The strip b n lies within the characteristic 
aft-cone from E n . Therefore, cr n and a n ‘ are the parts of S n lying 
respectively to the right and to the left of the wing. 

Let us return to the fundamental formula for the velocity potential, 
equation ( 4 . 2 ), which is in the characteristic coordinates 


cpp(x,y,z,t) = 


B ( x+ y) 

- J^e P 2 
2rt 




ScfXj_ 

e 2 

cos[a \|(x - i)(y - t]) - z 2 ] 

z) 

dz 

z=0 

\j(x - l)(y - tj) - z 2 


(6.24) 


In order to compute the velocity potential by means of this formula 
in those parts of the space (or, in particular, on the wing surface) for 
which the region of integration S(x,y,z) intersects the region S n of 


the x,y-plane, we must first determine 


6cp 

I 

6z 


t (*+ y ) 


outside the 


wing in the strips °i» °2’ 


• y G n , and. y cTp y • • • y y • • 


• y 


respectively. 
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Let us denote e ^ in the cr^, o^, • • • > cr n , . . . 

strips by 0, e< 2 >, e<3), . . 8 (“>. . . and in cx^ 1 , $2 > • • •> 

V> ... by 0', 0 f ( 2 ), . . ., 0'( n ), .... 

Let us construct the integral equation for 0^ . 

Let us express the velocity potential at the point N(x,y,0) in cr 2 
by formula (6.24) which is equal to zero everywhere in the strips a^ } 

°2> - • • °n (correspondingly in , a 2 ' , . . . <y '). 

Let us divide the region of integration into the three parts 
Sss + a+a-^'* as shown in figure 10. 


dcp-L -P 

The function e 

dz 


(x+y) 


2 = A(x,y) is given in s(x,y) on the 

^ p ( x+ y) 

wing. In c 1 , *(x,y) of a-\' , the function — ^ e 2 _ e'(x,y) is 

dz 

determined by the solution of equation ( 6 . 23 )- 


“P //-x\ 

In cr(x,y) we denote e 2 by 0^ ' (x,y) . Then we arrive at 

dz 

the integral equation satisfied by ©( 2 ) 


, (x+y) 



9 (2) (6,1) cos t>^ |)(y - n)] ^ ds 
/(x - I) (y - n) 


F^ 2 Ux,y) ( 6 . 25 ) 


26 


NACA TM 1383 


= r| = y and the known function F^^ is defined as 


where the limits of integration are hounded by Xg = £ = x and 

^ •*"» ^ ir n v%/^ 4*Vio Irvi ^T.TV> m s~f\ — I v(2) 


f< 2 > (*,y) = - IT Ad,,) cos N (x - O . fr . r _J t>l dr, di - 

. JU - £)(y - n) 

s(x,y) 

rr e'(M) cos N (x - sKy an as (6 . 26) 

J'J l/(x - EHv - Tl') 


t* 


We look for the solution of integral equation ( 6 . 25 ) in the form of 
the power series 


0^2) (x,y) = YZ. e ^ 2) (x,y) ^ 2n 
n=0 2n 


( 6 . 27 ) 


Moreover, by reasoning similarly to the preceding section we arrive 

( 2 ) 

at an integral equation for the coefficient 02 n in the general term 
of series (6.27) 



xzMt) 


4 n } (M) 


dn d| 


/(x - I) (y - ti) 


= 4 n^ x ^y) 


(6.28) 


where 


Fn 2 \x,y) = F n (x,y) + YZ f n^ k ( x >y) ( 6 . 29 ) 

k=0 

where, in its turn, 

f n 2)k(X ^ } = fan k)l 1 F (I ^ } [ (X " 0(y " l|) ] n “ k "2 *1 « 

V* (6.30) 

Equation (6.28) differs from equation ( 6 . 12 ) only in the form of the 
fC 2 ) function on the right side. Taking into account the condition on 
the boundary F n 2 ^ ( x E>y) = 0 for any n=0, 1, 2, . . . the solution of 

(6.28) for 9 p^ is obtained by using the solution (6.21) or (6.22) of 
( 6 . 12 ) as a final formula if F n ( 2 ) replaces F n in the latter. The 
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function Fn( 2 ^(x,y) depends on the coefficient 02k ^ where k=0, 

1, 2, . . ., n-1. Therefore, just as in the previous section, if the 
9 (2) for k = 0,1,2, . . ., n-1 are already found, then F n (2) in the 

right side of (6.28) is a known quantity. Therefore, the functions 
0 q( 2 \ 6 2 ^ 2 ^ • • •> e 2n^> * * * may he found successively. 

(2) (2) 

Let us note that F n ' , and therefore the coefficient 02n > 

depends only on the first n + 1 coefficients 0 O ', 02* > • • e 2n' 

of the series expansion of 


0'(x,y) 




e 


-pfei 

p 2 


in . 


eW, 

9 .( 5 ) 


Reasoning in the same manner, 
. . ., 0^. . . in 03, 


y 


e ,(M 


y 



we may find the values of 0'^, 

04, . . ., 0 N , . . . (correspondingly 

. in 0p*, 0 2 > • • • % )• 


Therefore, the velocity potential can he computed hy equation (6.24) 
at every point M(x,y,z) of the space for which the region S(x,y, z) 
intersects any number of strips 0jj or a^' . 


All the results hold for the case when the wing tips are not given 
hy one equation y = t(x) hut consist of curves given hy the equations 
y = i|f^(x) k = 1, 2, . . .,m . The same observation applies to the 

leading edges E‘E (or E-jE-^ 1 ) of the wing. Therefore, in our problem 

the wing contour may be piecewise smooth. 

If the frequency of oscillation a> of the wing he put equal to zero 

then the coefficients e Q , e Q ( 2 ), • • • > e 0^ N ^ . . • coincide with the 

values of the derivatives 8cp 0 /8z in the strips 0p, d 2 , • • •> . 

respectively, for the steady motion of a wing when the streamline condi- 
tion (1.6) on the wing is given in the form 


= A-j_(x,y) 
oz 
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7. INFLUENCE OF THE VORTEX SYSTEM FROM THE WING FOR A HARMONICALLY 

OSCILLATING WING 

1 . Let us consider the case when the region of integration S(x,y,z) 
in formula (k.2) for the velocity potential intersects the vortex sheet 
7 ^ as shown in figure 26(a) (see also fig. 11) . That is, let us consider 

the case when the trailing edge of the wing - the arc DI of the wing 
contour - or, so to speak, the vortex sheet, acts on the point M(x,y,z) 
of space. 

Using condition (1.10) we determine dq^/dz in the region ft of • 
the x,y-plane and shown in figure 11. 

The region ft is off the wing within the characteristic aft-cone 
from D and outside the characteristic cones from T. Therefore, ft 
is affected by the vortices trailing from the edge DT of the wing but 
not from D'T'. The region ft partially intersects the vortex sheet 
Si- 


Let us return to the characteristic coordinates x-j_, y^_, z-j_ which 

we introduced earlier by formula (6.3). 

As before, for brevity we omit the subscript 1 from the independent 
variables . 

Condition (1.10) fulfilled on 2 in the characteristic coordinates 
is 

+ u !a= + u = 0 (7.1) 

8t 8x 8y 

From equation (7.1) it follows that the function 

i<o x+y 

= c P i ( x >Y>°> t ) e u 2 

remains constant everywhere on the vortex sheet along lines parallel to 
the direction of the incoming stream, i.e., along vortex lines from the 
wing. 
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Since the velocity potential cpp = 0 everywhere in the x,y-plane 

off the wing surface and the vortex sheet, then it may he verified that 
cp^ possesses the specified property everywhere in ft. 

Let us construct the equation for the function 


Let us express cp^ at the arbitrary point N(x,y,0) lying in ft 
by using the basic formula for the velocity potential (6.2k) . We divide 
the region of integration S into three parts, as shown in figure 12, 
into s(x,y), c-^*(x,y) and cr(x,y) . The regions s and cr-j* are 

parts of the wing surface and defined above, respectively, which 

fall within the characteristic fore-cone from N(x,y,0). The region a is 
the part of ft in the same cone. The variables of integration in cr 
vary between xp <s | % x and X(|) % q y where xp is the abscissa 

of D and y = X(x) is the equation of the arc ET of the wing contour. 
The expression obtained for cp^ is differentiated in a direction parallel 

to the velocity vector of the impinging stream. 

Therefore we arrive at the integro -differential equation which d 
satisfies in ft 


*(x,y) = ~ e 

l? z Jz=0 



in ft. 







2 


dq d| = o(x,y) (7.2) 
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where (j. = -i 


2 2 

H — - a -- and. the known function is 
UO) 


®(x,y) = |- i - JJ A($,Ti)K^(S,n$x,y;X) dq d| - 
s(x,y) 

JJ 9(l,'n)K 1 (|,n;x,y;A) dTi d£ - 

°l(x,y) 

HA 2 JJ AU,T))K 1 (|,n;x,y;X) dri d£ - 

s(x,y) 

uA 2 JJ 6(6,Tj)K 1 (|,n5x,y;X) dn d£ 

di(x,y) 


(7.3) 


where K. U,tlj x,y;X) = CQS ^/( X , ~ ^ ( y and the operator 

fZ - 0 (y - n) 

iL. = The function 9 is determined from equation (6. 7) of the 

8L dx 8y 

preceding section. 

2. We will look for a solution of equation (7*2) in the form of the 
power series 


k2n 


(7 A) 


-3(x,y;A) = YZ ^ 

n =0 

Keeping in mind the absolute convergence of equation (7 .k) and using 
the expansion (6.8) for the cosine we obtain 


*(M;X)cos[X(/(x - |)(y - tj)] = 

00 k=n / , \n - k 


(-D 1 


n^H r , n, 

n=0 k=0 |2(n - k)J I 


[( x " 0 (y - t )] 11 


-k 


(7-5) 
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Substituting equation ( 7 , 5 ) , (6.8), and ( 6 . 9 ) into equation (7.2), 
the latter becomes 


00 k=n ^ i) 11 "^ 


f ffr x2n n r r 1 b* - «<* - d] 

^ So E5 l><“ - *)]* 


n-k- i 


dq d£ + 


k— n / s n-k 


A JJ t E WM) [<* - »<* - ^ 1 ^ « ♦ 


^ to fco [2( n - k )] 


»* Jf fl A2(n+1) F r !~ 1} , WM)[(* - 6)(y - T l)] n_k " 2 dq dl 

„ n=0 k=0 L 2 ( n " k )i 1 


8 

8l 



(-D n+1 

(2n)i 


A 211 


A(l,q)[(x - I) (y - q)] n 


1 

2 dq d| + 


. nn 00 l&=n n-k+1 

8 ff J— A 2n r~ (-1) 


8l 


a i 


* n=0 


£ rJ7 — -qv 9 2k (5 ’'' ) [ (x - s >k - ’>>] 

k=0 |_ 2 ( n " k )J 1 


n-k-| 


dq d£ + 


* (f E * 2(n+1) A(S,n)[(x - l)(y - 4" * d„ d| + 

J y 4— 1 (2n) • 


s ru=0 


k=n , . n-k+1 


n-k- | 


cr^ n=0 


E a2( “ +1) E k)] . <WM) [<* - » (y - D] - d, dl 


(7.6) 


Taking into account the uniform convergence of the series with 
respect to £ and q in both sides of equation (7*6), we integrate it 
term by term. Then, keeping in mind, the uniform convergence of the 
obtained series with respect to x and y which is also maintained after 
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differentiation, we differentiate the specified series term "by term with 
respect to x and y. After these operations on both sides of the 
obtained equation we equate coefficients in identical powers of A. There- 
fore we arrive at the integro-differential equation which the coefficients 
of equation (7»*0 satisfy 


f r r wm> _ dn di - - 

° x J Xd J x (£) /(x - £)(y - n) 


tr 1* wm) .. diai 

3y % J X ( i) Ax - |)(y - T|) 


= ®n(x,y) (7.7) 


where 


n+l 


®n(*,y) = ^rr- J- i Mm) 


dr] dS 


(2n) ! 8L 


I've - i)(y - t) 


, n+l 


^ ^ 7 vJJ A(£,ti)[(x - 0(y - n)] ^ dTj di + 


[ 2 (n - k)J 




[2(n - k)] . 
k=0 L J 


k=n-l ^ i) n-k 


E rail- 1 !) ! I “ 1 T %k(6 ' n) [ (x - 5)(y ■ ” ) l n ’ k ’ 5 4,1 45 + 


k=0 


[2(n - k)] 


h— n-1 . n— k+1 pp 

E /T wm>&. 


n-k- i 


k=0 


[2(n - k)]l 8L JJ 


l)(y - n)] 


dti d| 


(7.8) 
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in which the last sum and also the terms in p are defined for n>0. 

Let us note that the right side, $ n , of equation (7*7) f° r ^2n 
contains terms with coefficients but only for k = 0, 1, 2, . . . , 

n-1. 

Let us transform equation (7.7). We integrate by parts with respect 
to I the first integral on the left side of equation ( 7 . 7 ); "the second 
by parts with respect to tj, afterward we differentiate with respect to 
the parameters x and y, respectively. Equation (7*7) becomes 


X r y + 


2nt| ' 


XD x<?) (/(* - l)(y - n) 


dq d| = ^ n *(x,y) (7*9) 


where 


( x >y) = - 


1 py W x d^) 


dq + 


- X D Jx(x D ) l/y - TJ 
r -jiHil - l}d| + * n (x,y> 

Xd v< 


(x - 6) y - x(0 


( 7 - 10 ) 

Let us note that the first term in' equation (7*10) of the right side 
of equation ( 7 . 9 ) becomes infinite for x = Xp. 


Let us return to expression ( 7 . 8 ) for 4> n and separate out of it 

the terms corresponding to the value k = n in the first sum - the compo- 
nent 



0 2n ( ^> 


|/(x - D(y - T l) 


dn d| = R 


We integrate this integral by parts with respect to | keeping in 
mind that the limits of integration in a-j_* are x-g ^ I - x E and 
\jr(|) % Tj = y and that 0 2n (x E ,y) = 0. Then we differentiate with respect 

to x 
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R = 


y dn + 


C E 


/x - x d ^(x d ) /y - n 


X D /x - | 


r y 

l*U) ly . ,, 


an 


d| 


(7.H) 


Let us subject the desired function ■S in equation (7.2) to a sup- 
plementary condition. 

Let us assume that at the trailing edges - the arc DI (or D'T 1 , 
respectively) of the wing contour - and on the straight line DD* 

(figs. 11 and 12) - the intersection of the characteristic aft-cone from D 
with the z=0 plane (correspondingly the line D'Dp*) - the velocity of 
the perturbed flow, and therefore the function is a continuous func- 
tion, then the conditions are fulfilled 


d[x,x(x)] = a[x,'X(x)] 

= 9 L x D>y] 


(7.12) 

(7.13) 


These conditions are analogous to the Joukowsky condition for flow around 
a wing by an incompressible fluid. From equation (7.13) follows 


1 r v 1 °2n( x D> T1 ) / 

1 J - u dT] = 1 / u di] (7.14) 


l/x - x D J x (x D ) |/y - T) 


/x - x D ^t(x D ) /y - tj 


since x(xp) = i(xp) . 

Substituting equations ( 7 .H) and (7.14) in equation ( 7 .IO), the 
latter becomes 


tiW f«^WU n.x 

Jx D /(x - i)[y - XU] l d| 


(7.15) 


d| + 




X J) £ 


_1 

"77 


y ^(M) 


't(D 


dr) 


d.5 + V (x,y) 
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where 


V = " R 


(7.16) 


For n = Q, the right side in equation (7*9) is a known function of 
x and y 



^afeMO] 

fdX(S) n 

0 4 

/(* - l)'[y 

- x(s)]l d5 

P X E 

1 8 


8 o(i.i) d 1 

4 

/x - | dl 

y+u) 

/y - *1 


d£ 


d| - 


I // A(5,n) d ^ - 1: f 60 ( 1 , 1 .) 4,1 45 


|/(x - 0(y - t) 


/(x - £)(y - n) 

( 7 . 17 ) 


Let us solve equation (7-9) for ^2nx + ^2ny* 

The two-dimensional integral equation (7*9) is equivalent to two 
homogeneous equations 




d£ = ®n*(x,y) 


(7.18) 


and 


Jx “ \/F T 0 

r y a 2n p(M) + ^PrmCs^) , 

— dr) = d 2n *U,y) (7.19) 


'MO (7 - n 


each of which reduces to an Abel equation. Using the Abel inversion 
formula we find the solutions of equations (7.18) and (7.19) as 
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WW) = r 


1 ^n* (x D^ y) + l r V* (l ' y) 


31 fa - x d 11 Xd /x - l 


d| 


( 7 . 20 ) 


and 


«2n|(^y) + fl 2mr(S,y) = k 


2ny' 


! V*[S,x(0] ^ ^ ry *an n *(6jl) ai) 

J \ 


+ — 
X 


51 fa - x(i) 31 J *U) fa - n 


(7.21) 


Substituting equation (7.20) into equation (7.2l), first replacing 
in the latter by x, we obtain the solution of equation (7*9) as 


, $ n*[V X ^ x) ] 

«2nx( x ^) + ^2ny( x ^) = 3 7 = = 

|/x - x D l/y - X (x) 


1_ 1 

1,2 hr - x(x) 


- x $ n |*[^ x(x) ] 


X D 


d£ + 


_ L _i ^ y( x p^) 

3t 2 - x D -a (x) /y - TJ 




r y »n*T,*U^ 

( x ) /(x - l)(y - Tl) 


dq d| 


( 7 . 22 ) 


Integrating equation (7.22) along the straight line parallel to the 
free-stream between the limits of N(x,y,0) and N(x,y,0) we find the 
formula determining i3 in the general form of equation (7.4) 
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^2n( x >y) = <2n( x #y) + ■% 




2 1 

11 ^ x \/ x i - x d |/ x i + y - x - *( x i) 


dxi + 


' x f x l ® ri |*[^ x (xi)]d5 dx-L 


2 , — , , 
it ^x. 


: D ^xi -l^xi+y-x - x(x x ) 


_1_ 

2 

rt 


px^+y-x ^n , n*( x D> T l) dt l 


X(2t l J (/*! - X D / 


x-n |/xn + y - x - tj 


1_ 

2 

it 




px 1+ y-x ^nSt^ 6 ' 1 ^ 

x ( x i) \J x i - t ^ x i + y - x - t) 


dT] d| dx^ 

(7.23) 


If in equation (7.23) the coordinates x and y are taken as solu- 
tions of y-x + x- y = 0 and y - X(x) = 0 and the value of ^2n( x >y) 
is determined from condition (7. 12) on the trailing edge, then we find 
i3 p n on the vortex sheet. 


If in the same formula, the coordinates x and y are set equal to 
x = x D and y = y - x + x D and the value of •fl 2 n( x ^) is determined 

from equation (7.13) on the line x = xp, then we find •a 2 n outside the 

vortex sheet in the region it affects. 

Thus, through equation (7.23), we can compute successively the coef- 
ficients $ 2 , . . . , $2n> .... 


Therefore, the solution of equation (7*2) is constructed as the 
absolutely convergent series (7. 4) for any value of A. 

The coefficients •fl 2n ' axe expanded in the series 


•S’(x,y;A) = ^ ■d^ (x,y)A 2n 
n=0 


(7.24) 
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The function 


acp, 

•6' = e ^ in !)' (fig* 11) may be computed through 

dz 


equation (7*23) if the function X 2 ( x ) replaces x(x) in it (where 
y = X 2 (x) i- s equation of D'T' of the wing contour) and we inter- 
change the role of the coordinates. 


3* Let us consider the general case of the flow over an oscillating 
wing by a supersonic stream. Let the characteristic aft-cones from Eq 

and Eq' and Dq and Dq' intersect the wing as shown in figure 13. 

Then Eq (correspondingly E^ 1 ) , as shown above, are defined so that to 

the left on the leading edge equation (3*7) is satisfied and to the right 
it is not. The points Dq and Dq ' are, respectively, the most right 

and left points on the wing plan form. 


The space of the considered wing plan form as transformed by equa- 
tion (5*4) i s illustrated in figure l4. 

Let us divide the x,y-plane where the medium is perturbed into a 
series of regions : the regions considered in the preceding section, 

So Sq, . • • , S n , . . • , Sjj and the regions Aq, A£, • • • , 

Z!^, .... The region Sjj is the M-shaped region bounded downstream 

by the intersection of the characteristic cones from Dq and Dq' with 

the z = 0 plane. In the z = 0 plane, these lines are the upper 
bounds of the region of influence of the trailing vortex sheet. 

The region is M-shaped lying between the characteristic cones 

from Dn, D n ' , D n+ q, D n+ q 1 . We divide, in its turn, the part of the 

x,y-plane lying to the right and left of the wing, respectively, into 

the strips Oq, 02 , • • • cr n , . . . , defined above and into 6q, 

& 2 , . . ., 5 n , • • • and into aq', a 2 ', • • • > n n ', . . • , Ojj' 

defined above and 6q', 82', * • *, 5 n ', . . . correspondingly. The 

strip 8 n is that part of 2^ to the right and & n ' is the corresponding 

part of to the left of the wing. It is easy to see that the region ft 

defined at the beginning of this section is in 5q. 


In order to solve completely the problem of the flow over the wing 
shown in figures 13 and 14, the derivative c)cpq/c)z must be determined 
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in Bq, 62, 




and in 6 


I 

1 > 


5 2 ', 


J n > 


Let us denote the function 



and a', a'^, . . 

and 61 1 , 82 ' > . • • 


dcp 


1 -3 

e 


(x+y) 


dz 

^‘(n) 

• y v y 


-'ll y 


by *6^5)^ . . 

. in the bj_, 62 y • • 
strips , respectively • 


y 


Applying equation (6.24) for the velocity potential we construct cp^ 
for any point N(x,y,0) in 62* 


We divide the region of integration S which depends on the form of 

the function e 2 into the following: S = s + 0* + + S* + a, 

dz 

as shown in figure 15. This function is given in s. It was determined 
in cr* and a** in the preceding section by the solutions of equa- 
tions (6.7), (6.23), (6.27), etc. In s* it is determined by the solu- 


tion of equation (7-24). 


aq>l 

We denote e 2 in a by •a 

dz 


( 2 ) 


Using the boundary conditions (1.10) and (1.12) we arrive at the integro- 

( o) 

differential equation which a ' ' satisfies and which differs from equa- 
tion (7.2) only in the form of the right side. On the one hand the right 

side depends on the solutions 9, 0 . . ., 0^, 0', 0'^, . . • , 

0'(N) and on the other hand on the solutions a*. We construct a(^) i n 
the form of a power series in the parameter A. 


Requiring the fulfillment of equations (7*12) and (7*13) for a ^ 

we obtain for the coefficients aQ^, a 2 ^ , • . ., ^2sx'^ > • * * 

an expansion in series of a^^ of equations of the form (7 .9) which 
differ from each other in the form of the right side. 


( 2 ) 

The right side in the equation for the coefficient a 2n v ' in the 
general term of the series for a^) depends on the first n+1 coef- 
ficients of the expansion of 0^) and 9' where i takes all values 

( 2) 

less than or equal to N, and on the first n coefficients ’ , 

$ 2 . . ., a 2k (2) (te=0, 1, 2, . . ., n-l) of the series expansion of 


ho 
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the desired function Therefore, it is possible to find succes- 
sively the coefficients ■d 2 • • ., *2x1^ using the solu- 

tion (7.22) of (7.9) as a final formula if there is put in the latter, 
instead of 0 n *, right sides in the equations of the form of (7.9) for 

the respective coefficients of the expansion of fl Pr| (2) t 

By the same reasoning, values may he found of 
^ ^ • • . in • • • • 

Therefore the velocity potential may he computed hy equation (6.2A) 
at any point of the space perturbed hy the motion of the wing shown in 
figures 13 and 14 . In particular, the velocity potential may he eval- 
uated at any point of the wing surface. 

All the results are valid when the contour of the wing is piecewise 
smooth . 

If the frequency of the oscillations of the wing, <x>, he put equal 

to zero, then the coefficients %, • • •> ^0^, • • • coin- 
cide, respectively, with the values of dcp^/dz in &p, & 2 , . . ., 

5 k , . . . for steady motion when the streamline condition (l.6) is given 
on the wing as dcp^/dz = Ap(x,y) . 

We apply the proposed method of determining ckpp/8z for the oscil- 
lating motion of a wing hy constructing an integral equation, to wings 
of completely arbitrary plan form. For example, the wing contour may not 
he cambered hut may have the shape shown in figures 18, 2 k, etc. 

In all cases, the part of the x,y-plane where dcp^/dz must he deter- 
mined should he divided into the corresponding characteristic regions. 

Then successively passing downstream from one region to another, construct 
the integral and integro-differential equations using the boundary condi- 
tions on the x,y-plane. The solution of these equations for ckp^y^z or 

for functions related to 8cp^/8z is obtained as a series in even powers 

of the parameter A, which defines the frequency of oscillation. The 
whole problem of determining the coefficients of the expansion reduces to 
a double integral equation in each characteristic region. Each of the 
equations after transformation appears to he an equation of the same type 
which is solved hy means of a double application of the inversion formula 
for the Abel integral equation. The form of the wing contour is the limits 
of integration. The influence on the considered region, of determining 


NACA TM 1383 


41 


the desired function in the preceding upstream characteristic region, is 
reflected in the form of the function in the right side of the integral 
equations . 


8. FLOW AROUND AN OSCILLATING WING OF NON-ZERO THICKNESS 

1. Let us consider the motion of a thin wing at a small angle of 
attack (fig. 15 a) • 

Let the wing be moving forward in a straight line with the constant 
supersonic velocity u. Let an additional small oscillating motion be 
superposed on the basic motion of the wing so that the wing surface may 
be deformed. 

The normal velocity component on the upper surface of the wing will 
be considered given by 


°nn = Aou(*,y) + R - P * A 2u (x,y)e 1 ^ t ( 8 .l) 

and on the lower surface by 

°n l = A 0l ( x >y) + R,P * ^ (8.2) 

where A Qu and A Ql define the wing surfaces and 

A^ = A lu (x,y)e iot u( x ->y) and A 2 j = A 1 ^(x,y)e i(X l ( x >y) define the ampli- 
tude and initial phases of the additional oscillating motion of the wing. 
We consider the functions Aq u , Ap u and given at each point of the 

upper surface and Aqj, A-^, and given on the lower surface. The 

x,y,z coordinates were defined in section 1. 

The velocity potential cpp is 

cp p (x,y,z,t) = cp(x,y,z,t) + cp g (x,y,z,t) 


( 8 . 3 ) 
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The potential cp is specified hy the motion of an oscillating wing 
of zero thickness, which creates at each moment an antisymmetric flow with 
respect to the x,y-plane (fig. 15b). The potential cp B is specified by 

the motion of a thin oscillating wing with a profile symmetric relative 
to the x,y-plane. Therefore the motion proceeds in such a manner that at 
each moment the wing surface will be symmetric relative to a designated 
plane (fig. 15c). Such a wing creates a symmetric flow and cp s satisfies 


cp s (x,y,-z,t) = cp s (x,y,z,t) 


(8.4) 


Each of the potentials cp and cp s is represented, in its turn, by 


cp = cp 0 + ^ 


(8.5) 


<P S = lbs + 'fls 


( 8 . 6 ) 


where cp Q and cp 0g correspond to the steady motion of the wing and cp^ 
and cpp s correspond to the additional motion of the wing. 

Let us set up the streamline condition using the representation (8.3) 
for the velocity potential. 

We transfer the boundary conditions on the wing surface parallel to 
the Oz axis onto the projection £ of the wing on the x,y-plane 
(fig. 1). 

Therefore, we obtain the streamline conditions based on equa- 
tions (8.1) and (8.2) 
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8 z 


z=+0 


= AQu( x >y) + A2 u (x,y)e^ a * 


( 8 . 7 ) 


and 



8 z 


. z=-0 


= AQ Z (x,y) + R.P. 


AaiU.yJe 1 "* 


( 8 . 8 ) 


which must he satisfied on the upper and lower sides of z, respec- 
tively . 

Using equations (8.5) and (8.6) we establish boundary conditions 
for the desired potentials cpo> 9 l> TOs* and ^ls* 


Keeping in mind that on the z =0 plane the normal derivatives of 
the potentials cp^ s and cp]_ s are specified by the symmetry of the flow 

over the wing satisfying the condition 


Sep, 


Os 


Sz 




z=+0 


Sep, 


Os 


Sz 


z=-0 


Sep. 


Is 


Sz 


J z=+0 




8 z 


(8.9) 


z=-0 


We find the boundary conditions for ep^ s and ep^ s which must be satis- 
fied on the upper surface £ to be 

= R.P. r 2 (x,y)e ia3fc ( 8 . 10 ) 

z=+0 


.’^08 

= r n (x,y), ' 


8z 

(J 7 7 

rr L-H 

8z 


where the functions Fq and F^ related to quantities given on the 

wing surface through 

Aq (x,y) - A (x,y) A 2 u (x,y) - A (x,y) 

r 0 (x,y) = — r 2( x >y) 5 


(8.11) 


The conditions to be satisfied by ep Qs and cp^ s on the lower surface of 
£ are 


Sep 


Os 


Sz 


= - r 0 (x,y). 


Jz=-0 


Sz 


= - R.P. r 2 (x,y)e la3fc ( 8 . 12 ) 


z=-0 
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Since the normal derivative of the potentials cpQ and cp^ specified 
by the antisymmetric flow over the wing, on the z=0 plane, satisfy 


1 

& 


89 0 

- * 

8cp x ' 

► = . 

8cm 

dz 

z=+0 

dz 

7 

z=-0 

dz 

z=+0 

dz 


(8.13) 


the boundary conditions which must be satisfied simultaneously on the 
upper and lower surfaces of £ are 


89. 


dz 


2 = *o(x,y) 






3 R.P. k^_{ x.,y)el°& 


(8.14) 


where Aq and are related to quantities given on the wing through 


Aq = 


Aqu + ^ 1 


Ag = 


A 2 u + A2 1 


(8.15) 


The boundary problems for cfy(x,y, z,t) and cpQ(x,y, z) were set up 

in section 1 where in the case of a harmonically oscillating wing, equa- 
tion (8.14) rather than equation (1.6) should be taken on the wing. The 
solution of these boundary problems is contained in the present work. 

Let us formulate the boundary problems for cp^ s and 9 qs : 

I. Find cf>j_ s (x,y,z,t) satisfying equation (l.4), condition (l.ll) on 
the disturbance wave, condition (8. 10) on the plane region £ and 


8z 


0 


(8.16) 


everywhere in the x,y-plane off £ where the medium is perturbed. 

II. Find the function 93 S (x,y,z) satisfying equation (1.5), condi- 
tion (l.ll) on the disturbance wave, condition (8.10) in the plane 
region £ , and 


8z 


= 0 


(8.17) 


everywhere off z in the x,y-plane where the medium is perturbed. 
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Since the potentials q>L S and cp^ s are functions which are symmetric 
relative to the x,y-plane, it is sufficient to solve the problem for the 
upper half- space. 

The solution of boundary problem I is given by equation (4.2) . By 
means of this formula it is possible to compute the velocity poten- 
tial cp^ s everywhere since in the case of symmetric flow over a wing the 

derivative dcp^/dz is a given quantity for any point M(x,y,z) of the 
space in the region of integration S(x,y,z). To compute q>]_ s at M 
according to equation (4.2) the function 



dz 


z=+0 


must be substituted for 


*4 

Sz Jz=0 


R.P. r 2 (x,y)e ia5fc 


and integration is over that part of 


the wing within the characteristic cone from M. 

The solution of boundary problem II as is known (refs. 21 and 22), is 


8cp n 

given by formula (3.10) if the function <j 


is replaced by 


z=0 


SfOs 

t^J 


= r 0 (x,y) and integration is also over the region defined imme- 


z=+0 
diately above. 


If the wing is vibrating as a rigid body then the functions A2u and 
A 22 coincide and therefore, to solve the flow problem in this case, it is 

sufficient in antisymmetric streams excited by the motion of an oscillating 
wing with profile of zero thickness to superpose steady symmetric streams. 
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PART II 5 


To apply the integral equations method, explained in Part I of the 
present work, let us consider the problem of the flow over thin wings of 
finite span in steady supersonic flow. 


The velocity potential cp^, specified by the steady motion of the 
wing may be computed through equation (3-10) at those points M(xp,yi,zp) 
of the space for which the region of integration sCx-^y^ z-J , already 
known from Part I, does not extend outside the limits of the wing where 

8T0 . 

is given. 

dzp 

If 6tpo/8z2 appears to be unknown at any part of S, then, to use 

equation (3. 10) in these cases, where it has in the characteristic 
coordinates (6.3) the form 


c Po( x i>yi> z i) = " ~ 


2n 


s ( x i>yi> z i) 


'^0 

8zn 


dTi 1 di x 


, 1= 0 \j ( X 1 - &i) (yi - \) - Z 1 S 


(21.1) 


and to obtain the effective solution of the problem, it is necessary, first 
of all, to find dp 0 /dz 1 everywhere in S by constructing and solving 

an integral equation. 

1. INFLUENCE OF THE TIP EFFECT FOR STEADY WHIG MOTION 


1. The integral equation (5.1) in the coordinates (6.3) is, for the 
steady wing motion 


&i( x i,yi) 


drp-L dl-L 


d(x 1 ,y 1 ) 


/( x l - 61) (yi - 4 l) 


F( Xl , yi ) (21.2) 


5nie results of Part II, sections 1, 2, and 3 were completed in 
April, 1948 at the Math. Inst, of the Acad, of Science, USSR. 
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where 0-[_ is the value of 8cpQ/8zp on z , (fig. 6) and where the 
known function is 


F ( x l^l) = - 




drii dgp 


; ( X l' y l) 


/( X 1 - 61) (yi - nj 


(21.3) 


The function A given on the wing is 


Afa , yi ) . ^ - ug ° (it<y) - a 

' ^ F 8z n ^ k 


x i + yi . y i - x i 

2 ' 2k 


■} 


(21.4) 


dcpft 

It is easy to see that the velocity of the perturbed flow — — 

8z 

normal to the x,y-plane is related to ckf^/dzp through 


k ^b 


dz 


dzn 


The regions of integration in a are Xpg < |p < Xp and 
t(£p) < 4i < yp where, as before, y-j_ = t(xp) is the equation of the 
wing tip ED in the transformed coordinates and x n ^ is the abscissa 
of E in the same coordinates. The regions of integration for £p in 
s are the same limits ^ < ll < x x and ^(Ip) < m< ♦(fil) where 
y^ = ip^l) is 'the equation of the leading edge e'E of the wing contour. 


Let us note that equation (21.2) may also be obtained from equa- 
tion (6.5) if the frequency of the wing oscillation is set equal to 
zero in it. 


Let us delete the index "l" from the independent variables. 

We solve the double integral equation (21.2) with respect to &, by 

means of a repeated application of the inversion formula for Abel's inte- 
gral equation. 
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We write equation (21.2) as 


a 


x E 


l/x - I 


ry e (I,n) f tU) A(| n) 

't(S) /y - 4 *1^ /y - tj 


d| = 0 

(21.5) 


This is an Abel equation with right side identically zero, therefore, 
the brace equals zero for 5 = x. Hence, equation (21. 5 ) is equivalent to 


ry 0 (x,q) pi|f(x) 

— dq = - / dT , 

+ (x) /y - q ^tl(x) |/y _ Tj 


(21.6) 


which is also an Abel equation. Noting that the right side of equa- 
tion '(21.6) is, generally speaking, different from zero for y = ljr (x) 
we find the solution using the well-known inversion formula for the 
Abel equation 


9l(x,y) = i 


* /y - ■'l'(x) 


ri(x) 

4]_(x) 

A(x,q) 

fy - n 


r t(x) 


dq 


y = *(x) + 


* °r(x) /y - q dq 


A (- X ;V ) dq’ 

^l(x) ^ _ qt 


dq (21.7) 


Let us note that the solution (21. 7 ) for the steady motion of a wing may 
be obtained from the solution (6.22) of equation (6.12) for the vibrating 
wing if the index n and the frequency of oscillation 00 are both set 
equal to zero. 

Carrying out the operations specified on the right side of equa- 
tion (21.7) we find the solution of equation (21.2) to be 


e L (x,y) = - 1 


51 /y - t(x) J +i(x) 


^ (x) Jllr(x) — 

A(x,q) — L 1 dq (21.8) 


y - n 
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In a similar manner, we find the value = 0i'(x,y) in L * 1 

bz J 


(fig. 6) 


01 *(x,y) = - i 


1 1 


t 2 (y)/ x p 2 (y) - 5 
A(i,y) ' * 


* - ? 2 (y) , ^ 1 (y) /x - | 


d| (21.9) 


The functions x = y^(y) and x = ^(y) are, respectively, the equa- 
tions of the arcs ED and E*D* of the wing contour solved for x. The 

solutions (21.8) and (21.9) show that the velocity of the perturbed stream, 
when the axes ED and E'D' are approached from off the wing, goes to 

infinity as R ~2 where R is the distance of N(x,y, 0 ) from the points 

ED or E'D* (see fig. 7). 

2 . Let us find the velocity potential according to equation ( 21 . 1 ) at 
the point M(x,y,z) of space for which the region of integration S 
intersects the wing surface E and the region E3 or E?;'. 

The region of integration S in equation ( 21 . 1 ) is divided into three 
parts: S = s 1 + s 2 + s Q , as shown in figure l6 


cp Q (x,y,z) 



s 0 +s 2 


dil d| 

/(x - D(y - n) 


2 x 


ff &l(M) 

~ s i 


dn 

/(x - l)(y - n) 


(21.10) 


The limits of region s^ are x E < I = X A and - 1 l r (i)=n*y : 

./here x A is the coordinate of the point A which is the intersection 
of the characteristic forecone from M with the side edge ED of the 
..ring. The equation tj = y _ z 2 jx - l is the equation of the hyperbola 


I 


in which the aforementioned cone intersects the z =0 plane. The limits 
of region s 2 are x E ^ 5 = x A and tiU) - T - tU). 
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Using equation (21.8), let us evaluate the integral over s^_ in. 
equation (21.10) 


i = [ e^M) 


dn d(- 


/(x - |)(y - t) 


i rv y_ ^J 

rt(t) 

a U,V) \/iU) - V ln i 

dn d* 

* 'LJm) 

■Aid) 

Li 

l/n - t(l)(Tl - T)‘) 

i/(x - £)(y - n)-z 2 


( 21 . 11 ) 


we interchange the order of integration of tj 1 > *1 


i = - 


if /*%.,,■> 

* Jxe <A|tjU) /z - | 


f y - ti 

\Ji\ - i(l)0i - - 


z2 
x - E, 


- 1 


dt) 1 d| (21.12) 


The result of the inner integration is 

.2 


I* = 


r 

4(0 


X-l 


dt) 


<y- 


+(6)(n - u ’ )i / y - ~~ 2 ^ ~ ’i 


(21.13) 


I /tuT-n* yyT 


z 2 

z - V 


X - I 


Putting the value of equation (21.13) into equation (21.12) we obtain 


I = - // A(S,V) 


dn di 


/(x - i)(y - n) - z 2 


(21.14) 
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Equating (21.14) and (21.11) we obtain 


= -ff ACM) - - d1d - S — 

S 1 l^ x - 6)(y - t) - 62 /(x - l)(y - fj) - z 2 


(21.15) 

Therefore, to find the velocity potential, on the basis of equa- 
tion (21.1), at a point M(x,y,z) projected onto M' (x,y,0) in the 
x,y- plane as shown in figure l6, it is sufficient to integrate over Sq 


cp 0 (x,y,z) = - i // A(6,ri) ■ -? 1 ^ ■ ■■ - (21. l6) 

** £ l^x - !)(y - t,) - 


The limits of region sq are t]_(l) = q = y - z 2 /x-| and 

x A = ^ = Xg where is the abscissa of the point of intersection of 

the Mach forecone from M with the leading edge E’E. 

The velocity potential on the wing surface can be calculated from 
equation (21. l6) by setting z=0 in it and considering the region of 
integration to be x A Is £ a x and tq(£) - T - y because the lines of 

intersection of the characteristic forecone from M with the x,y-plane, 
in this case, are the lines i = x and tj = y. 

In order to compute the velocity potential at points of space, or in 
particular, on the surface of the wing for which the region of inte- 
gration S intersects simultaneously and that is, at points 

of space where there is felt the effect of both side edges ED and E'D‘, 
it is sufficient to integrate equation (21.1) over the region 
8 = S © + S ©, cross-ha 4ched region in figure 17. Hence the integral 

over S Q in equation (21.1) must be taken with the opposite sign, i.e., 
the plus sign. 
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3. Let us consider the wing of more general form shown in figure 18. 
Let the forward part of the wing have the "break, the arc EGG ' E-]_ 1 , in the 

wing contour which affects the flow just as do the side edges. 

Let us show how to compute the velocity potential at all points 
M(x,y,z) of the space disturbed by the motion of the wing, which is not 
affected by the trailing vortex sheet, in particular, on all points of the 
wing surface. 

We divide the wing surface into the characteristic regions shown in 
figure l8. 

If the region of integration S in equation (21.1) intersects 
regions 2, 2', 3 and does not intersect 4. then the velocity potential 
may be evaluated by using equation (21.16) (see figs. 16 and 17). 

The simple result which is expressible by equation (21. 16) does not 
hold in the general case. 

If S intersects k on the wing, in the curvilinear triangle K'OpK, 
then according to equation (21.1) dcffy/bz must first of all be found in 
the triangle. 

Let us express, by equation (21.1), the velocity potential at any 
point of K'O^K as equal to zero everywhere outside the wing and the 

vortex sheet from the wing, hence in K'OpK. Therefore, we arrive at an 

integral equation of the form of (21.2) for the function 

9^*(x,y) = ckp^/Sz in K'O^K but with a more complicated known function. 

Applying the Abel inversion formula twice, we arrive at the solution 
in the following final form: 


\ l 1 r* A(x, T)) \j 4r( x ) -i\ 

01*(x,y) = - J . ,- T / dT l - 

* /y - i(x) y - n 


A(|,y)\yi|» 2 (y) 


/x - V 2 (y) 


d| 


x - I 


(21.17) 


where y = \|r(x) is the equation of EG, y = \|>p(x) is the equation of E 1 
x = ? 2 (y) of Ei'G' and x = t 12 (y) of E^E. 
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Substituting equations (21.17), (21.8), and (21.9) into equa- 
tion (21.1) we obtain the formula for the velocity potential at M which 
has the projection M' shown on figure 18, and for which the region S 
intersects 4 on the wing and, therefore, the region K'O-jK outside the 

wing, as 


qb( x >y> z ) = - JJ A (M) 

S*(x,y,z) 


dri d| 


I /(x - S)(y - n) - Z 2 


1 JJ A <M) ^-1 1 

[**(£) - t(«][(x - l)(y - n) - z 2 ] 

* 2 s x * /(x - S)(y - n) - z 2 / 

1 If A !M) - tan' 1 1 

[vU) - n]|(x - 5)[y - **(£)]- z 2 | 

[f*(n) - * 2 (n)][( x - 5)(y - n) - z 2 ] 

1(2 s 2 * l/( x - £)(y - ’i) - z 2 1 

[t 2 ( 1 >) - 

|(y - n)[x - **01)]- z 2 j- 


dn d| + 


d| dTi 


(21.18) 


where y = i|r*(x) and x = t*(y) are the equations of GG' of the wing 
contour in terms of x and y, respectively. 

The region S* is the part of the wing shown cross-hatched in fig- 
ure 18. The regions S-]_* and S 2 * are part of S* and are marked in 

the same figure by horizontal stripes. The regions Sj* and S 2 * are 

bounded downstream by lines parallel to the coordinate axes passing 
through G and G* . The points G and G* are respectively the points 
with the largest x and y coordinate on the arc EGG'E^' . 

By combining the results of equations (21.1) and (21.18) there is 
found in the form of integrals taken over the wing surface, an effective 
expression for the velocity potential at points of space for which S in 
equation (21.1) intersects 5 or 6 on the wing and therefore A K'O-jK 

and and zy off the wing. 
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2. FLOW OVER WINGS OF SMALL SPAN 


1. Let us assume that the characteristic cones from Eq and Eq' 

intersect the wing as shown in figure 19 . This occurs, for example, 
for small span wings. 

Let us divide the x,y-plane where the medium is disturbed into the 
regions Sq, S-^, . . ., S n , .... 

The region S n is an M-shaped region lying between the character- 
istic cones from E n and E^ (or in one of them) and E n +q 
and E n+ i'. In its turn, we divide the part of the x,y-plane to the 
right and left of the wing into the strips 0q, 0 2 , • • • , cr n , . . . 

and o 1 l , a 2 ', . . ., a n ', • • ., respectively. The strip 0 n lies 

between the after cones from E n and E n _j_q. Therefore, 0 n is that 
part of S n lying to the right of the wing. The coordinates of E 
and E' with their indices are shown in figure 19 • The strip a n ' is 

similarly defined. 

Let the leading edge Eq'Eq be given as in part I, section 6, by 
the equation y = Vq(x) and the side edges ^i^n+1 anc ^ ^l'^n+l by 
y = \|r(x) and y = ^(x), respectively, or as x = i|r(y) and. x = i^Cy) 
correspondingly . 

To compute the velocity potential at M according to equation (21.1) 
in that part of space (or, in particular, on the wing surface) the region 
of which intersects S n of the x,y-plane but not Sj^q, we must first 

of all determine dcp 0 /dz off the wing in 0-^, a 2 , , o n 

and also in 0 -^', erg* ’ . . ., <%'> • • • • 

We construct the integral equation for &Pq/8z in the arbitrary 
strip 0 ^. 

Let us express a velocity potential which is equal to zero every- 
where off the wing and outside the region of influence of the vortex 
system from the wing, at N of the 0 ^ strip (fig. 20) according to 

the fundamental formula (21.1) 
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S(x,y,0) 




8z 


drj ds 


z=0 /(x - 5)(y _ ri) 


= 0 


( 22 . 1 ) 


The limits of integration in S axe x-j_ < | < x anri y^ < tj < y. 

For convenience in later writing, we make S a rectangle, which is~pos- 
sible since the medium ahead of the wing is not disturbed and dcfojdz 

is zero. The region S is shown in figure 20 bounded by the lines LN, 


NL^, L^Oi and O^L. 


Let us denote 
®2 > ' • •> ®k'> • 


z by 0 -l, 0 2 , 


> , 9^, • * * and 0^* , 


. in the respective regions a^, a 2 , . . . , a k , 

. and a 1 ', b 2 ', . . ., o k , . . .. 


In conformance with this new notation we write equation (22.1) as 


X D 


py Mm) 

i/ x - l |^t(fi) /y - ^ 


di) + 


-*U) 


A(i,n) 


^2(1) ft - 71 


dT] + 


i ^r 2 pyi +1 e^M) 


E 

i=i u yi 


/y - t 


dTj + 


r* 2^) (|,n) 

^k-1 ^ /7TT 


dn 


d| = 0 


( 22 . 2 ) 
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Applying the Abel inversion formula twice to equation (22.2) we find 
Gjj. for k > 2 


6 k (x,y) = -i_i— - 

* /y - nx) 


A(x,T))^(x) - T) 


%( X ) 


y - n 


dri + 


i= r^ 2 l~y i+1 V(*,l)ft r (x ) ptgM e k _ 1 '(x,r 1 ) ( V(x) - r, 

13 y - 11 4,1 + \_i ~ 


■dT] 


(22.3) 


Correspondingly, for 0^* we obtain 


V(*,y)=- £ 


11 ^/x - + 2 (y) 


' t2(y) A( 6 ,y)ft 2 (y) ~ 


't(y) 


X - 6 


d| + 


i=k^2 r x 1+ i e i (^y)^ 2 (y) - 1 ^ X 9 k-l^> y) /V y) ~ 1 

H J x x - | ^ + 4 k ! x - * 

i=l x i K_J - 


d| 


(22.4) 
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where the terms in equations (22.3) and (22.4) containing the summations 
are defined only for k ^ 3. 

If 0 2 , . • •, 9k_]_ and therefore, 0-^, 0 2 ', . . ., e k-l’ 

axe alxeady defined in a-^ 1 , a 2 ', . . ., cr^-l* then we can compute 

0k in crk for any k by means of equation (22.3). 

The value of in o^_ and c^' is determined by solving 

equations (21.8) and (21.9). 

The value of 8cp 0 /8z in d 2 is found from equation (22.3) by 
putting k = 2: 


e 2 (x,y) = - - 


1 1 


/t(x) - TJ 


i/y - t(x) Jy 2 (x) 


A(x,T]) 


y - t 


dr) + 


1 1 


^2(x) r t 2 (n) / i(x) - n \/i^(nj - l 

a(i,tj) L 1 d| dn 


it 2 /y - t(x) J y 1 ^+ l (t]) (y - n)(x - i)/x - f 2 (T]) 


(22.5) 


We find 8cf^/8z in cr 2 * in the same way 


e 2 '(x,y) = - h 


%(y) 


* /x - t 2 (y) J i(y) 


2U; , » \My) - 5 

A(l,y) 2-S d| + 


x - £ 


1 1 


J£ 2 /x - t 2 (y) Jx ± ^(1) 


>(y) p+(l) \/t P (y) - 6 \/Kt) - ti 

A(S,n) -^-= V - - drj d| 


(x - !)(y - n)/y - ^(i) 


(22.6) 


Thus, step by step we compute ckp^/dz in o k . 
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Using the solution of equation (22.3), we now prove the relation 


ft = 




dtj as 


8z 


z=° /(x - S)(y - T]) - Z 2 


= 0 


(22.7) 

where x-^* and Xg* exe any numbers satisfying x-^ < x^* = x A (^- A 
the coordinate of the point A shown in fig. 21), x ± < Xj* < x A . 


For the proof, we write 2 in the equivalent form 


ft = 




e^M) dT i 

. .. + 

/(x - l)(y - n) - z 2 




A(i,n) an at 

/(x - !)(y - n) - z 2 


i=k-2 X2 * y i+1 ei «(6,n)dndi ^ 

Jxj* J y ± /(x - l)(y - T l) - ^ 



e^’CM) d1 i d ^ 


'yk-i /(x - s) (y - n) - z2 


(22.8) 
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where 0 k in the first of the integrals is replaced by its value 
according to equatior^ (22.3). 

Then, we obtain 


a = - i f* 2 * [ Mt) 'jIlM - t. I. d | . 

* 9 ^, 2 (|) /x - l 


, r *2* r y± + i 0 ± '(M'VtU) - V 

-TZ I* an' d£ 

rt i=l Jxj* J y ± /x - £ 


p r * 2 * r * 2(6) e k _i(M')/+U) - n* 


71 ^x-l* 9y k _ x 


/x - I 


I* dll' d£ + 


r*2.* ri(t) A(M) dt] dg f 

x l* ^2 ^ ^ /(x - |)(y - t]) - z 2 


i= j^ r*2* /~' y i+l Qi'd^n) dii d£ 


i=l X l* y i /(x - g)(y - ti) - z 2 


p*2* r i2(l) ek-l'U^) dt! d£ 

J„ * J„_ . r. : 


x 2* y k-l /(x - |)(y - T]) - z 2 


(22.9) 


where I* denotes the integral (21.13) evaluated before. It is easy to 
see that all the terms in the right side of equation (22.9) cancel in 
pairs. Hence, equation (22.7) is proved. 

It is also clear that the following holds 


r\ V * n X 

/ 2 r y_Ti 

■^0 

' yp* '^x 1 

S Z 


df, dt} 


z=0 |/(x - 0(y - T}) - z 2 


= 0 
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where y^* and y 2 * are any numbers satisfying y-j_ = y^* < y-g and 
y l < y 2* “ y B ^ y B is tiie coor(iina ' te of B shown in fig. 21) . 

Using equations (22.3) and (22.4) it is possible to prove equa- 
tions (22.11) and (22.12) correspondingly 


<2* r y * 


an a§ 


x i* lSzJ z=0 /(x - 6 )(y - n) - z 2 


. x |*bl 

2 r m 1 5zJz=o 

31 ^ x l* ^ y l /(x - i)(y - n) - z 2 


tan- 1 /[( x - *)(y - nL - , z2 ^y* - t(o] di dn 
7 [(x - s)(y - y*) - z2 X+(s) - n] 


(22.11) 


where y* may depend on £ and satisfies \|r(x-i*) < y* = y — : 

- 1 - X-5 


r>y 2 * r.x* 

yi* ^¥ 2 (n) 


^ 1 dn d| 

^ z Jz=0 /(x - |)(y - tj) - z 2 


- 

_ g f y 2* r+atn) [~SI 


z=0 


y i Xl /(x - i) (y - Tj) - z 2 


tan' 


-1 /[^ x " ^ y " ^ " z2 ][ x * - * 2 (n)] 


a? an 


[(x - x*) (y - n) - z 2 ][? 2 ( t i) - 5 ] 

(22.12) 


where x* may depend on T) and satisfies ^(y^*) < x * = x 


y-n 


The relations (22.10) and (22.12) may be obtained, respectively 
from equations (22.7) and (22.11) if the role of the coordinates is inter- 
changed in the latter. 
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Let us note that the result of a single application of the 
Abel inversion formula to equation (22.2) or directly to equation (22.1) 
yields 



(22.13) 


Interchanging the role of the coordinates in equation (22.13) we obtain 



dj 

0 / x “ £ 


= 0 


(22.14) 


It is possible to consider equations (22.13) and (22.14) as rela- 
tions fulfilled along the characteristic lines LN and L r N‘ in the 
x y-plane where y and x are, respectively, the coordinates of N or 
jji lying off the wing and off the region of influence of the trailing 
vortex system (fig. 20). The points N and N' lie to the right and 
left of the wing, respectively. These relations can be useful for compu- 
tations . 

2. Let us turn to the fundamental formula (21.1). Using equations (22. 
(22.10) (22.11), and (22.12) we obtain, by calculation, the formula 

for the velocity potential at M(x,y,z) for which S intersects 

S n for any n > 0 


1 rr A(£,ti)an d! 

+ 1 IT 

AU,n)dr 

' 45 - 1 IT 

A(6,i 

dn di - 

_i 

jr 


|)(y - n) - z 2 

2* JJ JC 

s© * x 

- 0(y - 

n)-* 2 

/(x - 0(y 

- n) - z 2 

« 2 

H* 

±P } r* +i 

^2( X A) 


dT) d| 

x k=n-3 ,y k+1 

r+ 2 (yB) 

VCMXh. 

an 

d*- 

rt2 fc=l -X k 

•4(0 V(x . 

0(y - 1 ) - 

z 2 

' i 2 W- M 

A 2 (n) / (x 

- s)(y - n) - z 2 


1 

r* 2 W 

e n . 2 (S,n)a 2 


_i_ r*P( y B)] 

r+(y b) 

•a-2 <5 > n)n i 

at 

dn 

? k. 

Mi) Kx 

- 0(y - n) - 

z 2 

* 2 Jy a -2 

fix 

- 0(y - i) - z 2 




Adding 


: dr\ d| - 


( 22 . 15 ) 
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where the functions fi^ and fig are defined as 


ftl 


fi2 



/[(x - £)(y - n) - z 2 ] 

[y B - *(£)] 

[(x - l)(y - y B ) - z 2 

][iU) - tj] 

/[(x - S)(y - n) - z 2 ] [ XA - ^2^] 


[< 


x - x a) (y - n) 


z2 J [^2^) " *] 


and where the regions and are regions of the wing marked on 

figure 21. The region is the vertically- striped region on the wing 

surface. The region S 2 * is the horizontally- striped region of the wing 
surface. It is clear that S]_* and S 2 * intersect each other and S@ 
on the wing. 

The region lies off the wing and is vertically- striped in 

figure 21. This region is the sum of the regions over which are taken the 
integrals containing 9^’ for k=l, 2 , ... , n-2 in equation (22.15). 

The region S2 lies off the wing and is horizontally- striped in 
the figure. All the integrals are evaluated over it which together con- 
tain 0^ for ks=l, 2 , . . ., n-2. 

If M is such that S in the basic formula intersects S n falling 
in the characteristic cone from En and lying outside the cone from Eh', 
then n must be replaced by n-1 in the second sum and in the last term 

of equation (22.15). If S falls inside the cone from E n * and lies 

outside the cone from En then n-1 must be substituted for n in the 
first sum and the penultimate term of equation (22.15). 

Let us note that the sums in equation (22.15) are defined for n > 3 

and the last two terms in equation (22.15) for n > 3. 
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If n=l, then the formula for the velocity potential in equa- 
tion (22.15) is limited to the first two terms. This result was already 
obtained before. 

If n=2, the formula in equation (22.15) is limited to the first 
four terms, the region of integration is shown in figure 22 . 

Thus, to evaluate the velocity potential, by equation (22.15), at 
a point M(x,y,z) which has the projection M’(x,y,0) shown in fig- 
ure 21, it is necessary, first of all, to compute 0^ for k=l, 2, 3> 

. . ., n-2 by equation (22-3) for k>2 and by equation (21.8) for 
k=l (0^' correspondingly). 

As an example we present the expression for the potential for n=3 
in the expanded form 


) = - A. jj A(6 <’»> d ' 1 d t. + JL ij A(i,il)dj_d| 


s» A* - l)(y - n) - z 2 


S0 


/(x - 5)(y - 1 ) - z 2 




1 dtj d§ - A. I — 1 ' - \— l/J “ 2 


A(g,n)fl 0 dn as 


n2 S-L* /S - l)(y - Tj) - Z 2 * 2 Sg* /(x - i) (y - TJ) - z 2 


_1 P^2 ^ X A ) r+Ol) n*U) A(|,T)»y7U) - 


Jy 2 o x . 


1 * l(i) /n - Hi) (n - v)/(x - t)(y - n) - z 2 


drj’d^ dt) + 


^(yu) rt 2 (5) r t 2 (rj) 


ACi^n) 


*3 


\l 


^ (n) - t' ^ 


yi 


* l(,l) h - * a <ti) U - i')\/(x - 5 )(y - 


dg'dTj d£ 


ri) - z 2 


( 22 . 16 ) 
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The region of integration in the last two integrals over | and r] 
are, respectively, the regions Sj_ and S£ lying off the wing and shown 

striped in figure 23. 

Formula (22.15) for the velocity potential contains an n-iterated 
integral with the integrand an arbitrary given function on the wing: 
dcpo/d z = A(x,y) . 

In the general case, it is not possible to reduce the number of 
iterations in the computation of equation (22.15) for arbitrary wing- 
tip shapes since the arbitrary functions \|r, ^2, and A all contain 

the variables of integration. If the f met ions \|/ and are fixed 

then the wing to be considered has completely determined tips and it 
is easy to see that all the integrals in equation (22.15) are reduced 
to double integrals taken over the wing surface with integrands containing 
the arbitrary given function A(x,y) which defines the form of the wing 
surface . 

Let us turn to the wing of small span which has a break in its 
leading edge as shown, for example, in figure 2k. 

The derivative cfcpo/6z may be evaluated in cr^ and by equa- 

tions (21.8) and (22.3)- It is impossible to evaluate depo/dz in cr^ 

using equation (22.3) and, therefore, a surface -integral equation must 
again be constructed which will also reduce to two Abel equations but 
with more complex right sides than occurred for o ^ in figure 19* 

Hence, we note that it is impossible to construct one formula which 
would determine dcpo/dz for all cases, but a single method of solution 

may be shown to depend on the wing plan form. 

The formation of the surface -integral equation for depo/dz is 

explained above, for each characteristic region. Each of these equa- 
tions is of the same type, reducing to two Abel equations with different 
right sides in different cases. In particular, the right side of one 
of the Abel equations, in some cases, may be identically zero. 
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3. INFLUENCE OF THE VORTEX SYSTEM FROM THE WING FOR STEADY WING NOTION 


1. To study the influence on the air flow of the trailing vortex 
system in steady motion, it is convenient to operate with the acceleration 
potential $q which, in linearized theory, is related to the velocity 

potential derivatives in the characteristic coordinates through 


\ 


'0 = u | c P0x + 


(25.1) 


Let us turn to the wing shown in figure 25. Let us take a point 
M(x,y,0) on the wing surface, which lies between the characteristic 
cones from D and D'. Therefore the trailing edge DT affects M. 

Using equation (21. 15) the velocity potential at M according to 
equation (21.1) is 


cRo(x,y,0) 


i_ CP A(s,n) an d| _ _i_ IT a(i,q) dg dg 

/(x - i)(y - n) 2jt 'jJ i/(x - 0(y - ti) 

s=S]_+Sq ' B 2 

( 23 . 2 ) 


where the regions s = S]_ + s Q and s 2 are shown in figure 25. The 
region s 2 belongs to G, considered in section 7 of part I and shown 
in figure 11. We denoted the derivative dcpo/dz in ft by -9 where 
this derivative is an unknown. 

We subject 8cpo/dz to an additional condition, analogous to the 

Kutta-Joukowsky incompressible- flow condition. We assume that the 
perturbation velocity potential at the trailing edge - the arcs ET 
and D'T' of the wing contour (figs. 11 or 25) - and therefore, the 
specified derivative, is a continuous function. Then the respective 
conditions are fulfilled: 
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fl[x, x(x)] = a[x,x(x)] ( 23 . 3 ) 

*[x,5fe(x)] =A[x,X 2 (x)] (23.4) 

where, as above, the function y = x(x) is the equation of DT and 
y = )^(x) is the equation of D’T‘ of the wing contour. 


In order to obtain the acceleration potential $q at M on the 

wing surface, we must take the derivative of equation ( 23 . 2 ) in a direc- 
tion parallel to the oncoming stream. Before differentiating the double 
integral with respect to x and y we integrate by parts - in the first 
case with respect to |, in the second with respect to q. 

During these operations, we use equation (23.3) and the relation 
( 22 . 13 ) which is fulfilled along characteristic lines, and which on the 
line DD* (fig. 25 ) is 



yv'i „ 

, — a, » 

I y - 1 


■X(x D ) A (: 


-D> 


/y~ 




dq 


( 23 . 5 ) 


We keep in mind, moreover, that the limits of integration of S]_ are 
X D = ^ = X A *(£) = T = ^ 1 ( 5 ) where xp is the abscissa of D and 

X A = x A(y) is the abscissa of A, the limits of s Q are x A < £ < x 
and ii(5) tj < y and finally the limits of s 2 are x D < | < x A and 

x(6) < n < y. 

After the specified operations, the results of differentiation are 


<fb x (x,y) + <iby( x >y) 


- i ff 4. ( M) + V s ’l dn a6 . 

2lt ° /(* - 6)(y - 4) 


s l +s 0 


2n 


s 2 


*|(M) + ^( 1 , 11 ) 

/(x - S)’(y - tj) 


dq d£ - (23.6) 


_1_ 

2rt 


p A [1,^(1)] 

J 

f, 

. -L 

1 \ 

f(x - i)| 

1 1 

1 

*l(s)l 

1 1 

i « j 




( 23 . 6 ) 
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where the arc l = HP is shown in figure 25 • In. order to evaluate the 
acceleration potential o 0 at M according to equation (23.6) it is 

first of all necessary to determine + -8 in s 2 . 


2. Let us construct the integral equation for d x + d y . Let us 

express the accleration potential through equation (21. l) at an arbitrary 
point N(x,y,0) outside the wing in a affected by the vortex sheet 
trailing from the wing 


^(x,y,o) = _ _L 

2jt 


A (M) 


(x,y) - S)(y - 1) 

<r(x,y) K* - t)(y- 1 ) 


dT] dS 


nr 


dq d? 


( 23 . 7 ) 


for which the limits of integration in a are x D < | < x and 

*U) = *1= y and in s, £ varies between the same limits but q 

between ^(g) < r\ < x(g) (fig. 26). 

Let us differentiate this expression in the free-stream direction. 
Since, according to the _ condition ((l. 10) of part I) the velocity poten- 
ia cpo off the wing in the x,y-plane remains constant along lines in 
the specified direction, then the left side of equation (23.7) goes to 
zero as a result of differentiation and therefore we obtain 


8x 


'X(I) 


A (M) 


/X D 


_8 

r 

8x 

x d 

8 

r 

dy 

J X D 

8_ 

r x 

dy 

X U 


/(x - |) (y - qj 

fl(M) 


dq d| ¥ 


x (£) /(x - 0(y - t) 
fX(5) A(|,q) 


dq dg + 


J 


^l(^) /(x - £)(y - q) 
^ ’ 3 (|, q) 


dq dg + 


X(l) S)(y -T) 


:dq d| = 0 


( 23 . 8 ) 
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We integrate the first two integrals in equation (23.8) by parts 
with respect to after which we differentiate with respect to x. The 
result is 


a_ r x f m a(m) 

8x j J*i(S) /(x - |)(y - Tj)" 



r x ( x c) a(x d ,t) 

^x(x) i/y - n 


an 




/y - n 


•dil 


d| 


and 


(23.9) 


_8_ 

8x 



/" y 8(M) 

/(x - |)(y - 11) 


rdq d£ = 


/ x - X D 



— — — — dTj + 

/y - n 



JX(|) 


d(M) 

n 


dn 


■d| 


(23.10) 


Keeping equation (23.5) in mind, which is fulfilled on the characteristic 
DD* we substitute equations (23. 10) and (23*9) into equation (23.8) 
obtaining 



1 



r dMidn.i r* l) ab&« + 

'-'xd) /y - n /y - n 




x(t) /y 


dT) + 


x(i) 


A(i,n) 


^iU) /y - n 


dn 


_d 

8y 


■ d| = 0 


(23.11) 
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This equatipn is equivalent to 


_8_ 

8x 


_8_ 

dy 




d(x,n) 

dn + A 

r x(x) a(x,,) 

/y - n 

Sx 

^l(x) /y - n 

d(x,Ti) 

dri + A 

r * x) a(x,ti) 

i/y - Tj 

dy 

^+i(x) /y - n 


(23.12) 


according to the inversion of the Abel integral equation. 

We integrate the last two integrals in equation (23.12) by parts 
with respect to t) after which, as above, we differentiate with respect 
to the parameter. Using equation (23*3) we arrive at 



* x (x,t|) + 3(x,n) 

— 1 an = - 

/y - n 


rx(x) A x (x,q) + A T] (x,n) 
fl(x) /y - r] 


d,r\ - 


{x,* x (xj] 

dii(x) 

1 

y - ^(x) 

dx 


(23.13) 


Let us apply once again Abel's inversion formula, keeping in mind 
that the right side of equation (23.13 ) , generally speaking, is different 
from zero for y = X (x) we obtain the solution for d x + dy as 


d x (x,y) + \( x ,y) 


1 1 

51 /y - x (x) 



- Ax(x,n) 




Atj(x,ti) 


|/X(x) - T1 

y - n 



1 

fy - x(x) 




d%(x) 

dx 


1 : /x(x) - i x (x) 
y - ^(x) 


(23.14) 
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Using equation (23. 1*0 we prove 


s 2 


+ * n U>q) 

^(x - 0(y - q) 


dq d| 



A (£,q) + A ( | , q ) 

» dq d| 

/(x - l)(y - q) 


r a[i,*(d] 

lii(l) 

1 - 

z i /(x - i)[y - tiU)l 

d£ 

k 


(23.15) 


where Z-j_ = RQ,. The regions s 2 and s-^ are shown in figure 25. 


Substituting equation (23.15) into equation (23.6) we obtain the 
formula for the acceleration potential 


$ 0 (x,y) 


pp Ag(6,q) + Aq(M) 
— = %* + -toy = - S JJ ^===- 




2rt 


/(x - £)[y - tqU)] 


1 - 


d^q(|) 


d| 


-d£ 


(23.16) 

where L = QP, the direction of the integration is shown by the arrows in 
figure 25. 

Thus to evaluate the acceleration potential at M on a wing stir- 
face two integrals, the surface integral over Sq and the contour inte- 
gral over L of the leading edge are to be computed. 

Let us turn to equation (23.12) and write it in the form 


_6_ 

r - 

8cp 0 

dn + _8_ 

8x 

-'tl(x) 

dz 

V. i 

i z =o /y - ^ 


r ■ 

^0 

dq 

^i(x) 

'dz 

J 

z=o /y - q 


= 0 


(23.17) 
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Interchanging the role of the coordinates in equation (23.17) we 


obtain 


8_ 

r - 

^Pol 

dt + a 

r . 


8x ■ 

^ ti(y) 

hz )z=0 

/ x - t 

^♦ x (y) 

dz 


A S = 0 

z=0 /* " 6 


( 23 . 18 ) 


where x = ^(y) is the equation of E*E of the wing leading edge solved 
for x in terms of y. 

It is possible to consider equations (23.17) and (23.18) as rela- 
tions which hold along characteristic lines in the x,y-plane where the vor- 
tex sheet has effect. 

Relation (23.17) is fulfilled along characteristic lines parallel 
to the Oy-axis (the line NQ, on figure 26) ; the y-parameter is the 
ordinate of a point lying off the wing to the right, in the effective 
range of the vortex sheet (point N in fig. 26). Relation (23.18) is 
fulfilled along lines parallel to the Ox-axis; the x parameter is the 
abscissa of a point lying off the wing to the left. 

If the point N is thus located to the right of the vortex line DH 
or to the left of D'H’, then along characteristic lines the respective 
relations (22.13) and (22.14) also hold. 

If N is located to the left of DH or to the right of D'H’, 
respectively, then relations (23.17) and (23.18) hold along characteristic 
lines. In this case, equations (22.13) and (22.14) are not fulfilled. 

In this section, we wrote down the transformation and obtained 
the formula for the accleration potential in the simplest case of the vor- 
tex sheet affecting the flow. 

For any other case, the potential <j>q is found in an analogous way. 
In each case an integral equation is constructed for + •a y . All the 

integral equations are of the same type but with different right sides in 
the different cases, and they are inverted by means of a double application 
of the Abel integral equation inversion formula. 

In the following paragraph we present results defining the accelera- 
tion potential <£q at any point of a wing surface. 
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3. Let us find the velocity potential cpo(x,y, z) at a point M 
lying within the characteristic aft-cone from D and outside the charac- 
teristic aft-cone from D* . The region of integration S in the funda- 
mental formula (21.1) intersects the plane region 0 , (fig. 11) in this 
case . 


The projection M 1 of M on the x,y-plane is shown in figure 26a. 

Starting from condition (1.12) (of part I) we express the derivative 
&Pq/c)z for any point where the velocity potential equals zero and where, 

simultaneously, the effect of the vortex sheet is felt through the same 
derivative at points located upstream on the same characteristic line 
with the point studied. To do this we reason just as we did to obtain 
formula (21.8). We then obtain the desired representation for the 
derivative 



Using equation (23. 19 ) it is easy to prove 


(23.19) 



y- 


fdqp, 


0 


1+y^lSzJz^ 


O 


at) dt 

/(x - 5)(y - n) - z 2 


px 2 * 

pl+y D -x D 

5 c p 0 ] 

dT) d| 

* 

X 1 

°VlU) 

az Jz =0 , 

/(x - l)(y - tj) - z2 


(23.20) 


by the same methods used in proving equation ( 21.15). 
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The limits of integration in equation (23.20), x^* and x 2 *, are 

any numbers satisfying Xp = Xj_* = Xp and Xp = x 2 * = Xp where Xp is 

the coordinate of the point F shown in figure 26a. The point F is the 
intersection of the vortex line DH, which has the equation 
y = x + yp - x D , with the characteristic cone from the point with the 

coordinates (x,y,z). 


In particular, there holds 
dn d| 


rr 

cxp 0 

JJ 1 

dz 


S 2 


fa® 


z =° /(* - I)(y - <1) - x 2 £ 


0 

- ^ • z=0 /(x - i) (y - T)) - z £ 


dT) d£ 


( 23 . 21 ) 

where the regions S-j_ and S 2 are shown in figure 26a. The region S-j_ 
is marked with horizontal and the region S 2 with vertical crosslines. 

Keeping in mind equation (23.21) we obtain an expression for the 
velocity potential at the point M defined above 


cfh( x ,y, z) = - — 
2rt 



A(6,n)dTi d£ 

/(x - I) (y --tj) - z 2 


— IT d(i,T 1 )dn di 

s' /(x - i)(y - n) - z 2 

(23.22) 


where Sq and S' are shown on figure 26a. 

Therefore, the region of integration S in equation (23.22) inter- 
sects the wing surface only in that part of 0 which lies to the left of 
the vortex line DH. 


Before evaluating the velocity potential by equation (23.22) it is 
necessary to determine dcpo/8z = d in the region S' of ft. 

We find d from the solution (23. 14) if the_latter is integrated 
in a free stream direction between N(x,y) and N(x,y) . Hence in order 
that the obtained expression correspond to the value of the deriva- 
tive 8 cPq/8z = i9 in 0 , to the left of DH, the coordinates x and y 
on the vortex sheet should be taken as the solution of the equa- 
tions y-5J-y D +xp = 0 and y = X(x) and the value of d(x,y) Is 

determined from equation (23.3) at the trailing edge. 
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If the x and y coordinates are set equal to x = x D and 
y = y-x+x D and the value of d(x,y) is determined on DH from the 

solution of equation (21.8) then the obtained expression will correspond 
to the value of 8cf^/8z in Q. to the right of DH off the vortex sheet 

but in its sphere of influence. 

4. PRESSURE DISTRIBUTION ON A WING SURFACE 

1. Let us consider a wing of arbitrary plan form. Let the wing 
contour in the characteristic coordinates be given by the following equa- 
tions: The leading edge E'E by y= i|f(x) or x = f^y), the side 

edges ED and E'D' by y = t(x) and y = t 2 U) or x = Y(y) and 
x = ^(y) , the trailing edges DT* and D'T* by y = x(x) and 
y = )<2(x) or x = X(y) and x = X^y) . 

Let us find the pressure of the flow on the wing surface. 

According to the Bernoulli integral, the pressure difference of the 
flow above and below the wing is related to the acceleration potential Oq 

by 

p(x,y) = P z (x,y) - p u (x,y) = 2p$> 0 (x,y) (2^.1) 

where p is the density of the undisturbed flow. 

We divide the wing surface into the ten characteristic regions shown 
in figures 27 and 28. 

Let us express the stream pressure on the wing surface in each 
characteristic region by the function A(x,y) which is given on the wing, 
defining the shape of the surface. 

We denote by M and M with a subscript the ends of line segments 
parallel to the coordinate axes and lying in the x,y-plane. It is clear 
that these segments are parts of the lines of intersection of the charac- 
teristic cones, with vertices in the x,y-plane, and the x,y-plane itself. 

Region I is the region where the tip effect is not felt. This part 
of the wing lies ahead of the characteristic aft-cones with vertices at 

E 1 and E . 
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Region II is where the tip effect is felt but not the influence of 
the trailing vortex sheet. This region lies between the characteristic 
aft-cones from. E' and E and D and D' . At M of region II, for 
which the lines and MgM^ intersect on the wing as shown on 

figure 27 , the pressure difference is 


t \ up 

P(x,y) = 

3t 


D(S,n;x,y)dT} dg 


Si 


—■ JJ~ D(g,il;x,y)dTi dg + 

S 2 


up 

it 





1 - 


L 


dg 



d-t(y) 

dy 


J B[i(y),r];x,y|dTi - 

L 1 


up 

Jt 


Jl 


di|r 2 (x) 

dx 


y' B [l,+ 1 (x);x,y]dg 

L 2 


( 24 . 2 ) 


where is the region of the wing bounded by the lines MM]_, MMg, 

and $>2 the region bounded by MpM^, and the 

arc L = and where 


Ag(g,n) + 

/(x - g)(y - tj) 


B(£,^;x,y) 


A(l,n) 

i K* - t)(y - n) 


D(g,ii;x,y) 
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If the lines and do not intersect on the wing, as shown 

in figure 28, then the pressure difference is 


p(x,y) = “ If JJ D (M; x >y) dT l d 6 - ^ / B [6,^(6) jx,y] Jl 


-%.ir 

D(i,n; 


s 2 


up _ 

it 

1 - 

4 ^(y) 

qy 

UP _ 

1 - 

dt 2 (x) 

or 

-J- *"■ 

dx 


d+i(l) 


d l 


d£ 


Ll 


B[+(y),Ti;x,y]dTi - 
B[|,t 2 (x);x,y]d| 


(24.3) 


L 2 

where is hounded hy the lines MMg, and 

L = M^. 


Arrows in the figures show the direction of integration in the con- 
tour integral and the integrals taken over the lines L-]_ = and 

L 2 = M^IV^. 

In region III, which lies between the characteristic cones from E 
an d the characteristic cones from E 1 , D and D', the pressure differ- 
ence is 


p(x,y) = - ^ II D(|,ti;x,y)dti d| - ^ / B [6,^(1) jx,yj -1 1 - 


up 


S 1 


dtpU) 


up 

3 t 


1 " } J B L*to'^ x >y] dTl 


as 

(24.4) 


d| - 


The pressure difference in region III* is expressed in the same way. 
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p(x,y) 



D(g,T^x,y)dn d£ 





(24.5) 


Region IV lies in the characteristic cones from E and E* and D 
and outside the characteristic cone from D' . Region IV* is defined cor- 
respondingly. At M(x,y) of region IV, when M]_M^ and intersect 

on the wing, the pressure difference is 


p(x,y) = - — II D(i,n;x,y)dT) d| + ^ II D(£,Tj;x,y)dTi d| + 


up 


Si 


s 2 


^ j B [1,^(1) ;x,y] 


L 



dti(|) ' 

~ 


J /^B[|,t 2 (x);x,yjd^ (24.6) 

L 2 

For the M, for which MpMj and M 2 M 4 do not intersect on the wing, 

the pressure difference is expressed by equation (24.5)* Similarly, the 
pressure difference for region IV 1 is 


up 

it 


— - 1 - 


df 2 (x) 
~dx 


p(x,y) = " -7/ / D (Mi x ,y) dT l d £ + 77 II B(^,n;x,y) dT i d£ + 


s l 


S 2 




2£<i _ 

si 


d t(y) 


< 3 y 


/ B 

L n 


[^(yM^y]*! 1 ! 


(24.7) 
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if and MyM^ intersect on the wing. If these lines do not inter- 

sect on the wing the pressure difference can he expressed by equa- 
tion (24.4) . 

In region V, which lies within the characteristic cones from E, E', 
D and D* where the influence of the trailing vortex sheet is felt, the 
pressure difference is 


p(x,y) = - JJ D(£,il;x,y)dT] d£ + ^ JJ D(|,r];x,y)dii d| + 


L 


S2 
dt X (l) 


di 


- d| 


(24.8) 


if M]_Mj and M 2 M 4 intersect on the wing, and 


p(x,y) = - — JJ D(i,iijx,y)dTi d| 
S 1 


b IMiUW] 

L 


1 - 


<%U) 


d| 


■ d| 


(24.9) 


if they do not intersect. 

In region VI, lying in the characteristic cones from E and D 
a pri outside the characteristic cones from E* and D* (also in 
region VI') the pressure difference is expressed by equation (24.9). The 
pressure difference for ra§ ion I has the same form. 

Thus, if M, at which the pressure is desired, is in one of the 
regions II, IV (IV 1 correspondingly), or V, as shown in the figures, then 
to set up the regions and contours of integration in the pressure formulas 
it is necessary to proceed as follows: Draw two lines MM-j_ and MM 2 

upstream from M to intersect with the side (or trailing) edges of the 
wing. From these points of intersection Mp and M 2 again draw lines 
Mq_Mj and NyM^ upstream to intersect the leading edge E'E at Mj and 

M4. 


If M is in region III or VI (ill 1 or VI' correspondingly) then from 
M draw the lines MM^ and MM-|_ upstream; the line MM^ immediately 

intersects the leading edge E'E at M^; MM 1 intersects the side edge 
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ED in the case of region III or the trailing edge DT' in the case of 
region VI. From the point of intersection Mp again drav the line 

to intersect the leading edge E'E. 

Let us consider particular cases. 

(i) Let the side edges of the wing ED and E'D* be straight 
lines parallel to the free stream. In this case 

_ 1 

8y dx 


and, therefore, formulas (24.2) and (24.3) are simplified substantially, 
because the last two terms in them become zero. 

A particular wing of this class is the rectangular wing. 

(il) Let the wing surface be such that 


D (M; x >y) = 0 


This holds, firstly, when the wing surface is a plane, i.e., the 

function A = u ^0/k is given on the wing, where 0o is the angle of 

attack, as a constant. 

Secondly, this holds when the wing surface is linear, generally 
speaking, uncambered, with generators lying in planes parallel to the 
y = x-plane (x, z-plane in the original coordinates), then the derivative 
of the function A(x,y) given on the wing satisfies the rela- 
tion A^. = - Ay. In particular this is a wing with a cylindrical surface 

formed in the manner described. 

In these cases, only the oontour integrals and the integrals over the 
line segments Lp and L2 remain in the formulas for the pressure. 

(ill) The pressure formulas take an especially simple form when 
the wing surface is such that the function D(&,rj;x,y) =0 on the wing, 
at the same time as the side edges ED and E'D' are straight lines 
parallel to the stream (combination of cases I and II). In this case, the 
pressure difference above and below the wing in any region can be repre- 
sented by 


P(x,y) = ± J Bft^D^y] 
L 


1 - 


<%(!) 
1T~ 


• d| (24.10) 
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where the plus sign is taken if the lines and Mr>M^ intersect on 

the wing and the minus sign if these lines do not intersect on the wing. 

Hence, the pressure on the wing surface is expressed by the curvi- 
linear integral taken over the arc L of the wing leading edge. 

(IV) Let the wing plan form be such that the points D and E 
and E* and D’ coincide. In this case, calculation of the pressure 
on the wing surface is also simplified because there are no regions II, 
III and III' on the wing. In particular, the trapezoidal wing belongs 
to this case. 

2. The pressure formulas show that there can exist a geometrical 
locus F*(x,y) = 0 where the pressure on the wing p(x,y) = 0. Down- 
stream of this geometrical locus, the pressure difference p = p^ - P u 

is negative. 

For example, if D(£,T)jx,y) = 0 on the wing then the geometrical 
locus F* = 0 is found in the region of the wing lying inside the char- 
acteristic cones with vertices E and E' and passing through either 
regions II and IV or through IV and V or or lying entirely in V. 

The first case occurs only when K, the intersection of the lines OpK 

and O 2 K parallel to the coordinate axes, appears to be outside the 

region of influence of the vortex sheet, as shown in figure 27 , for 
example. In all these cases, the points T and T 1 are on the geo- 
metrical locus of F* = 0. The curve F* = 0 may also be shaped convex 
downstream and not as shown on the figures. 

Let us write the equation for the geometrical locus F* = 0. 

In region II: 


nHv) 

d-fpU) 

d| 

° ^2 [^2 ^ x il 

1 d| 

V J 

/(x - !)[y - +i(l) ] 



1, ajfrll 4 - _ J 

d\|r 2 (x) 

/x - ^ 1 [i 2 ( x )] 

ay J F \/ x - t(y) 

dx 

r \l y - t 2 (x) 


(24.11) 
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In region IV: 


F*(x,y) 


rx(y) 

at-, U) 
1 - ^ 

as 

4i[t 2 (x)] 

d| 

/(x - S)[y - tpU) ] 


dilr 2 (x) J I x - t[t 2 (x) 
2 [ 1 dx J '/ ][y - t2( x ) 


(24.12) 


In region V: 


F*(x,y) = tpjXgU)] - x(y) = 0 (24.13) 


If the side edges of the wing are lines parallel to the free stream 
direction or the wing is such that E and D (E 1 and D’ correspond- 


ingly) coincide, then F* = 0 takes 
not changed, but in regions II and 
of equations (24.11) and (24.12) 

F* = ?i(> 2 (x)] 

and 

F* = tiftjjfr)] 


, simple form. In region V it is 
IV, we have, respectively, in place 


- t(y) = 0 

(24.14) 

O 

II 

lx 

1 

(24.15) 


In all cases when the pressure difference on the wing, according to 
equations (24.2) to (24.9), is expressed only by means of curvilinear 
integrals taken over L of the wing contour, it is easy to construct the 
zero-pressure curve graphically, keeping in mind that the zero-pressure 
curve in these cases is the geometrical locus of such points M on the 
wing surface for which the points M 3 and M 4 on the wing contour coin- 
cide. That is, the arc on the leading edge over which the curvilinear 
integral is taken shrinks to a point. 

We construct the zero-pressure curve as follows: From each point 

Nq on the leading edge we draw the lines NqN]_ ^ Nq^ parallel to 

the coordinate axes intersecting the side edges ED and E’D* as shown 
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in figures 29 and 30, or the trailing edges as shown in figures 31 and 32. 
From the points of intersection Np and N2 within the wing again we 

draw lines NpN* and N 2 N* parallel to the coordinate axes. The geo- 
metrical locus of N*, where these lines intersect, is the desired zero- 
pressure line. 

For example, for a symmetric wing, if the side edges ED and E'D’ 
are parallel to the stream, the zero-pressure curve passes through G 

and G' and is the line equidistant from the leading edge (fig. 30) . 

The points G and G’ are shown on figures 29 to 32. If E and D, 

E* and D', correspondingly, coincide and the trailing edges are straight 
lines then F* = 0 passes through G and G* and is the curve obtained 
by inverting the leading edge E'E relative to the center of inver- 
sion 0*. The center 0* is the point of intersection of the trailing 
edges (fig. 31). 

If the wing is asymmmetric and if the side edges ED and E'D* are 
parallel to the free stream then the zero-pressure curve is the reflection 
of the curve equidistant to the leading edge and passing through G and 
G’, relative to the line equidistant from the side edges (fig. 29). If 

the points E and D, and also E* and D’, coincide and the trailing 

edges are straight lines making identical angles with the stream then the 
geometrical locus F* = 0 is the reflection of the curve obtained by an 
inversion, with center 0*, of the leading edge and passing through the 
points G and G' relative to the line equidistant from the side edges 
(fig. 32). 

3. All the obtained results are generalized to the case when the 
leading edge E’E is given not by one equation y = ^(x) but consists of 

segments of smooth curves given by y = tik(x), where k = 1, 2, . . ., n 

with n any integer. In such cases the surface and contour integrals in 
the formulas for the pressure should be divided into component parts for 
the actual evaluations . 

The side, ED and E*D', and trailing, DT' and D*T, edges may 
also be piecewise smooth. 

The same generalization holds for the previous three sections. 

A. All the results are generalized in the case of the asymmetric 
flow over a wing which occurs, for example, in the motion of a yawed wing. 

Let us consider a wing of arbitrary plan form with an angle of 
yaw 7 as shown in figure 33- 
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The pressure on the wing can be computed by the same formulas if 
the equation of the arc Eq'Eq, in the coordinates transformed to the 

origin 0, is taken as the function y = ^(x) . 

The equation of EqDq (correspondingly Eq'D 0 ') is y = i(x) . In 
this case EqDq acts as the wing tip. 

Finally, for the trailing edge, DqTq, we have the equation y = X(x) 
( correspondingly for Dq 'Tq ' ) • 


5. As is known, knowing the acceleration potential or the velocity 
potential on the wing surface, we can easily compute the aerodynamic 
forces on the wing. 

In order, we represent the aerodynamic-force formulas using the ori- 
ginal coordinate system shown in figures 1 and 2. 

The lift P on the wing is 


P = 2p JJ <D 0 (x,y) dx dy 


(24.16) 


where the region of integration in Z is defined by ^o(y) = x = x i(y) 
and y D ' < y < y D where x = \|r 0 (y) is the equation of D’E'ED and 
x _ x-jiy) is the equation of the trailing edge D'TT'D (figs. 27 and 28) 
The limits yp and y^ are respectively the coordinates of D* and D 
of the wing. 

Since according to linearized theory 3>o( x > y ) = u then 

integrating (24. l6) over x and keeping in I^’ 

tial is zero on D'E’ED from conditions (1.11) and (1.12; oi pari 1, 

the lift is 


P = 2pu 




% 


[x L (y),y] 


dy 


y D-, 


If the trailing edge is piecewise smooth, then in actual computa- 
tions the contour integral must be divided into its component parts. 
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The egression for the moment Mq^ due to lift relative to the 
Oy-axis is 


The moments relative to the other axes have the same form. 

6. The explained theory can he generalized to the case of the flow 
over a tail or over a biplane in tandem. 

We proceed as follows to obtain formulas to compute the pressure on 
the tail taking into account the influence of the wing. 

Express cpox + TlQy at M(x,y) on the tail using the basic formula 
(21.1) • In the expression for Cty-) X + q^y under the integral sign inseri 
•3 X + on the vortex sheet. The function d x + -3 is found from the 
Abel integral equation which is constructed by the method of section 3 • 

In the case of flow over the tail the different characteristic 
regions on the tail must be separated just as was done in figures 27 and 
28 for the uniform motion over a wing. 

Only in this case, to divide the tail surface into regions, there 
must be taken into account, on the one hand, the wing effect and on the 
other hand, the tip effect and also the effect of the vortex sheet of 
the tail itself. 



dx. dy 


(24.17) 
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APPENDIX 

EXAMPLES 


The following examples, solved by N. S. Burrow and M. M. Priluk, 
will serve to illustrate the methods explained before. 


A. Arrow-Shaped (or Swallowtail) Wing 


Let us consider the arrow-shaped (or swallowtail) wing plan form 
where the leading edges are f ormed by the segments AD and AD ' and 
the trailing edges by the segments DB and D'B as shown in figure 34 . 
Let the following geometric parameters he given: the angle between 

the leading edge and the free -stream direction; 62 the angle between 
the trailing edge and the free-stream direction and l the wing semispan. 

The equations of the wing leading edges in the x,y characteristic 
coordinates with origin at 0 are 


line AD 

yi = 


1 + cot a* tan 5p 
line AD' 


<(1 - cot a* tan 5p)x^ + 21 cot a*l 


Vi = 


1 - cot a* tan 5^ 
and the trailing edge equations are 


J(1 + cot a* tan 5p)x^ - 21 cot a* 


line DB 


yi 


1 

1 + cot a* tan 



cot a* tan 5 2 ) x i + 21 cot 



line 

yi 


D'B 


1 

1 - cot a* tan 



+ 


cot a* tan 



21 cot 



where the angle a* is the semiapex angle of the characteristic cone. 

Let us consider the wing for which 6p > a* and 5 2 > a*; that is, 
a wing surface not affected by the trailing vortex sheet. 
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We will assume that the wing surface is a plane inclined by an 

&P n 

angle p 0 to the free-stream direction. Therefore, the derivative — ^ 

will be a constant everywhere on both sides of the wing surface and will 
be given in the form 


— = - u(3q tan a* (Al) 

oz 

In conformance with the method we divide the wing sin-face into the 
three characteristic regions la, lb, and Ic, with each region having its 
own analytic characteristic solution and taking into account the angular 
point A of the leading edge (fig. 3*0 . Let us compute the stream 
pressure on the wing surface in each region. 

Using the formula ( 5 . 9 ), we find the pressure in the regions la and lb, 
lying outside the characteristic cone from A, to be 


p = U ^— = 2 u 2 p 3 0 tan a* 

W 71 


(A2) 


This formula shows that the pressure in regions la and lb is a constant. 

In region Ic, lying inside the characteristic cone from A, we find, 
by using the same formula, the pressure to be 


p(x,y) = 


2u 2 pPq tan 5p 
^Jcot 2 a* tan 2 8p - 1 


1 tan' 

rt 


o , ll + cot a* tan 5p l cot 8p - xp 
— tan _J - 1 


1 - cot a* tan 5p l cot &p - Xp 
1 t cot a* tan y-^ - 1 cot 5^ 

(A3) 


[1 - cot a* tan 6p yp - 1 cot 8p 
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In the original coordinate system shown in figures 3 h and 35, (A 3 ) 
"becomes 

, , 2u 2 pPo t 911 8 1 

p(x,y) = 


Jcotfi a* tan2 5^ - 1 

f 


2 /l - cot a* tan 8]_ l cot 8p - x + y cot a* 

it ^a 31 y 1 + cot a* tan y cot a* + x - l cot 8^ 


o i /l + cot a* 

^ tan -K/ 

* V 1 - cot a* 


tan Sp 2 cot 8p - x + y cot a* 


tan 5]_ y cot a* + x - l cot 8^ 




These formulas show that the pressure is constant along each ray from A 
in region Ic. 

Shown in figures 36 and 37 , respectively, are the pressures along a 
section ApB-j_ parallel to the y-axis and along the section AgBg par- 
allel to the x-axis. 

The lift P of the considered wing is 


2u 2 pp 0 l 2 (tan 8 X - tan & 2 ) [ ( 2 + __;l / cot a* tan &p - 1 
P ' tan Scoots a* - l\ + * “ V cct a* tan 6l + 1 


2 tan & 2 _ - tan & 2 ^ /cot a* tan 8 p + 1 

it tan &p + tan & 2 V cot a* tan &p - 1 

tan 3 8 2 


tan 8 i(tan 2 8 p - tan 2 82 ) 


tan 


_2_ /cot a* tan & 2 - 1 


'cot a* tan 8 2 + 1 


(A5) 


The lift coefficient C z is 


C z = 


"t 811 8 l 


y|cot 2 a* tan 2 8 p - 1 


o n /cot a* tan S-i - 1 
1 - Z tan" 1 - \ / + 

it 


cot a* tan 8 p + 1 


2 tan 8 i ~ tan 5 2 ^ ^ — 1 jcot a* tan 8p + 11 
11 tan 8^ + tan & 2 Vcot a* tan &i - l[ 

l6 Pq tan 2 &2 


it (tan 8 p + tan 8 p)\/cot 2 a* tan 2 8 2 - 1 


tan’ 


_p cot a* tan & 2 - 1 
cot a* tan & 2 + 1 


(A6) 
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As is well known, the wave drag coefficient C x is related to the lift 
coefficient through C x = |3 qC z . 

Let us consider particular cases of (A6) . In the limit as 5 q_— > 
we obtain for the triangular wing 

C z = 4p 0 tan a* (A 7 ) 

the well known result^ for the lift coefficient of a triangle. 

Comparing (A6) and (A 7 ) we conclude that for identical wing speeds 
and identical angles of attack the lift coefficient of the arrow-shaped 
wing exceeds the lift coefficient of the triangular wing. 

In the particular case when 62 = 81 , we obtain the infinite span 
arrow-shaped wing. In the limit as & 2 — > 5]_ (A 6 ) yields 

4p Q tan 5 X 

\jcot2- a* tan2 - 1 

This result shows that the lift coefficient of an infinite span arrow- 
shaped wing equals the lift coefficient of an infinite span slipping 
wing with slip angle Sp. 

Formula (a 6) shows that with increasing 5-^ and &2> the angles 

between the leading and trailing edges and the free stream, respectively, 
the wing lift coefficient decreases. The dependence of C z for an 

arrow-shaped wing on and &2 is shown in figures 38 and 39 • 


B. Semielliptic Wing 

Let us consider the wing plan form which is half an ellipse as shown 
in figure 1+0. Let the semiaxis a-j_ and b^ of the ellipse be given. 

Let us assume that the wing moves, as shown in the figure, in the direc- 
tion of the axis of symmetry. 


^See the work of M. I. Gurevich: On the Lift of an Arrow-Shaped 

Wing in Supersonic Flow. Prik. Mate. Nekh., Vol. X, No. 1+, 1946. 
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The equation of the leading edge , the line D'D, in characteristic 
coordinates with origin at 0 is 


y l = “ *1 

and the trailing edge equation in these same coordinates is 
(a-[_ 2 - bj_ 2 cot 2 t 2a 1 h 1 cot a 


yi = 


a-^2 4- b-^2 cot2 a* 


In the original x,y coordinates the trailing edge equation is 


= + hhLi£ 


(Bl) 


The plus sign relates to the arc CD of the ellipse and the minus sign 
to the arc CD ' . 

Let us assume that the wing surface is a plane inclined at an angle 

8cp n 

Pq ‘to the free -stream direction, therefore the normal derivative — ^ 

d Z;L 

as given by (Al) . 


Let us consider the flow around the semiellipse when the character- 
istic cones from D and D' intersect on the wing surface. In con- 
formance with the method we divide the wing surface into the four 
regions I, VI, VI', and V. Region I is outside the characteristic cones 
from D and D ' , hence the vortex sheet trailing from the wing exerts 
no effect here. Region VI is within the characteristic cone from D 
but outside the cone from D'. Conversely, VI' is within the cone 
from D' and outside the cone from D. Region V, however, falls within 
both the characteristic cones from D and D ' . 


Using the formulas, we compute the pressure in each region on the 
wing surface. The pressure in I is constant everywhere and expressed 
by (A2). In VI the pressure distribution in the x,y coordinates is 
given by 


p = u 2 pPq tan a* x 


o -1 cot a* ILy + B f + 2a, b.. cot a*\/R. - f 2 

q _ sin 1 2 l 1 1 VI I 




xB^ 


(B2) 
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where 



f-j_ = x + y cot a* 



Similarly for region VI*. The pressure distribution in V is 


p(x,y) 



x 


sin' 



+ 



(B3) 


xB-^ 


where f 2 = x - y cot a* and B-]_, B 2 , and are as defined in (B2) . 

Graphs of the pressure distributions along the respective sections ApB]^ 
and A 2 B 2 parallel to the y-axis are given in figures 4l and k2 and 
along the corresponding segments A^B^ and A^B^ parallel to the 
x-axis are shown in figures 4-3 and kb. Spanwise section lines A^B-l 
and A^Bg are shown in figure 4-5; whereas chordwise section lines A^B^ 
and A^B^ are shown in figure 40. 

If the semiaxis of the ellipse are given in a special way; namely, 
if there exists between the semiaxes the relation a-]_ = b-^ cot a*, then 

formula (B2) for the pressure distribution in region VI simplifies, 
becoming 


Ukis corresponds to the case where the characteristic cones with apexes 
at D and D* intersect the wing trailing edge on the axis of symmetry 
of the wing at the point C; consequently the region V on the wing now 


p(x,y) = u 2 pp Q tana*<l - 2. sin-1 


cot a* y + ^2a^ 2 - (x + cot a*y) 2 


x 


(B4) 


vanishes. 
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In the general case for the flow around a semielliptical wing, it 
may "be shown that on the surface of the wing in region V, there exists 
a certain curve along which the pressure difference between the upper 
and lower surfaces of the wing reduces to zero. Downstream from this 
curve on the surface of the wing the pressure difference becomes nega- 
tive. We find the equation for this line of zero pressure by equating 
the right side of (B3) to zero. 


After obvious transformations, we represent the desired geometric locus 
in the following final form 




l+a-^^bp^ cot 2 a* ^ 


ap 2 - b^ 2 cot 2 a* 


) 2 ( 




(B5) 


where 






(B6) 
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0^00 * and h, . In order that the zero-pressure line should not pass 


Constructed in figure A 6 is an isometric view of the pressure on a 
semielliptic wing in the general case when (B7) is not fulfilled and 
there exist the regions I, VI, VI', V on the wing. 


Let us consider the wing of hexagonal plan form shown in figure 47. 
Let the leading edges he the lines OE^, and OE^', the side edges E^D 

and Ex'D* parallel to the free stream, and the trailing edges DB and 
D'B. In characteristic- coordinate space, the wing has plan form as shown 
in figure 48. 

Let us assign the following geometric parameters: 0 — the angle 
the leading edge makes with the free stream; 7 - the angle the trailing 
edge makes with the free stream; l - semispan and h chord. 

Let us consider that wing for which 0 > a*, 7 > a*. The first 
inequality means that the wing surface extends outside of the character- 
istic cone from 0. The second inequality means that the wing surface 
is outside the sphere of influence of the trailing vortex sheet. 




(B7) 


C. Hexagonal Wing 


NACA TM 1383 


93 


The equations of the lines forming the wing contours are: 
line OE-j_ 


y = x tan a 


or in characteristic coordinates 



where 


1 - cot a* tsLn a 
1 + ctg a* tg a 


here m < 0, since 5 > a*; the line 0E-j_* 


y = - x tan a and = mx-j^ 


the line E-^D 


y = l and yi = + 2 cot a * J 


the line E^_ 1 D 1 


y = - l and y^ = x^_ - 2 cot a*l 


the line DB 


y = - x tan 7 + h tan 7 


and 



+ n x 


the 
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and finally D 1 B 


where 


y = x tan 7 - h tan 7 and y-j_ = + rig 


1 + cot a* tan 7 

itn = ~ 

1 - cot a* tan 7 


2h cot a* tan 7 
n-, = 

1 + cot a* tan 7 


2h cot a* tan 7 
^ 1 - cot a* tan 7 


In conformance with the method we divide the wing surface into the 
13 characteristic regions shown in figure 48 . 

Assuming that the surface of the wing is a plane, we give the stream- 
line condition in the form (Al) and we compute the pressure in each 
characteristic region. We produce below the results of computing the 
pressure on the wing surface as formulas already transformed back to the 
original coordinate system. 

The pressure in la and lb is constant and expressed by (A2) . 

In Ic the pressure is 


p(x,y) = 


2u 2 pP q 


rt \] - m cot a* 


— + tan“l 
2 


-i 

m f x + 


- cot a* y 
cot a* y 


- tan 


-1 


- m 




- cot a* y 


* yl 


cot 


a* y| 


(Cl) 


Hence it follows that the pressure is constant along each ray starting 
from 0 in Ic. In Ilia 


p(x,y) = 


2u pp Q (l - m) 
it i /- m cot a* 


tan 


-1 


2 cot a* (l - y) 


(m - l)(x + cot a* y) + 2 l cot a* 


(C2) 
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In Illb 


p(x,y) = 


2u pP 0 (m - 1) 
Jt \l - m cot a* 


- tan 


-1 1 


- m 


/x - cot a* y 
j x + cot a* y 


tan 


-1 


2 cot g*( i - y) 


(l - m) (x + cot a* y) + 21 cot a* 


tan -1 - cotgy. 

f x + cot a* y J 


(C3) 


In IIIc 


p(x,y) = 


2u 2 pp 0 (l - m) 

jt j/- m cot a* J (l - m) (x + cot a * y) + 2mZ cot a* 


tan 


2m cot a*(y - z) 


(<*) 


In Ila 


p(x,y) = 


2upp 0 (l - m) 
it ij- m cot a* 


tan^ - 


/ 2m cot a*(y - z) 

(l - m)(x + cot a* y) + 2mZ cot a 


tan - ^ i /(l ~ m ) (x - cot a* y) - 2Z cot <x* 
I 2 cot a*(z + y) 




(C5) 
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In lib 


2u pp 0 (l - m) f 1 [ 

x - cot a* y 

P(X ' y) = , fTS c ' cT ^ T=J— cot a. y ‘ 


tan 


-1 


2 cot a *(z - y) 


(l - m) (x + cot a* y) - 2Z cot a* 


n /(l - m)(x + cot a * y) - 2Z_cota* 
^ an / 2 cot a*(Z + y) 


-1 1 — — /x - cot a* 

tan x f- m ,/ — - — 5 

' (X + cot a* 


(C6) 


In lie 


p(x,y) = 


2u PP 0 (m - 1) 


jt - m cot a* 


tan 


_ 1 1 /( 1 - m)(x - cot g* y) + 2mZ cot g* _ 


F 


m 


2 cot a*(z + y) 


tan 


-1 


1 x - cot a* y 


- m 


/ 


x + cot a* y 


+ tan - -*- m 


x - cot g* y 
x + cot a* y 


tan - -*- if- m 


2 cot a*( Z - y) 


(l - m) (x + cot a* y) + 2mZ cot a* 


(C7) 


Formulas for the pressure distribution on the wing surface in 
regions Ilia', Illb', IIIc', and Ila* may be obtained from (C2) , (C3), (C4), 
and (C5) , respectively, if coordinates appropriate to the specific 
regions are chosen. 

The formulas for the pressure show that there is a zero-pressure 
line on the wing surface, downstream of which the pressure difference 
below and above the wing becomes negative. This line is formed of the 
two segments KN and KN* the equations of which are 
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y = x tan 5 - 21 tan a* tan a y = - x tan 5 + 21 tan a* tan a (C 8 ) 

and which are parallel to the leading edges E^O and E^'O. 

The zero-pressure line may easily he constructed graphically. 

Graphical representations of the respective pressure distributions 
in the sections ApBp, A 2 B 2 , A 3 B 3 , A]^, and A 5 B 5 parallel to the 

y-axis are given in figures k 9 , 50, 51, 52, and 53. 

An isometric pressure surface is shown in figure 5 k for the 
hexagonal plane wing. 


Translated by Morris D. Friedman 
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Figure 11. 
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Figure 14. 
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